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historische Entwickelung der Methoden darzulegen. DarQber hinaus aber 
mars anf eine eingehende, mit Beweisen versehene Daratellimg, wie sie 

selbständigen, voa umfangreichen Quellenatudien unabhängigen Ein- 
dringen in die Disciplin erforderlich ist, auch hei den breiter angelegten 
Eeferaten der Deutschen Mathematiker -Vereinigung, in welcher haupt- 
sächlich das historische und teilweise auch das kritische Element znrj 
Geltung kommt, verzichtet werden. Eine solche ausföhrliche Darlegung, 
die sldt mehr in dem Charakter eines auf geschichtlichen und litterttr 
rischen Stadien gerundeten Lehrbuches bewegt und neben den rein 
wissenschaftlichen auch pädagogische Interessen berücksichtigt, erscheint 
aber bei der raschen Entwickelung imd dem Umfang deS zu einem 
grofsen Teil nur in Monographien niedergelegten Stoffes durchaus wichtig, 
zumal, im Yergleiche z. B. mit Frankreich, bei uns in Deutschland die 
mathematische Litteratur an Lehrblichem über sp(!zielle Gebiete der 
mathematischen Forschung nicht allzu reich ist. 

Die Verlagsbuchhandlung B. G. Teubner giebt sich der Hoffiiimg 
hin, dafs sich recht zahlreiche Mathematiker, Hiysiker und Astronomen, 
Geodäten und Techniker, sowohl des In- als des Auslandes, in deren 
Forschungsgebieten derartige Arbeiten erwünscht sind, zur Mitarbeiter- 
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sehafb an dem üntemeluneii entschliefsen möchten. Besonders nahe liegt 
. die Beteiligung den Herren Mitarbeitern an der Encyklopädie der Mathe- 
matischen Wissenschaften. Die imifangreichen litterarischen und speziell 
fachlichen Studien, welche für die Bearbeitung yon Abschnitten der 
Encyklopädie vorzunehmen waren, tonnten in dem notwendig eng be- 
grenzten Eahmen nicht vollständig niedergelegt werden. Hier aber, bei 
den Werken der . gegenwärtigen Sammlung, ist die Möglichkeit gegeben, 
den Stoff freier zu gestalten un^ die individuelle Aui^assung und Bichtung 
des einzelnen Bearbeiters in 'höherem Mafse zur Geltung zu bringen. 
Doch ist, wie gesagt-, jede Arbeit, die sich dem Plane der Sammlung 
einf&gen lafst, im gleichen Mafse willkommen. ' ^ 

Bisher haben die folgenden Gelehrten ihre geschätzte Mitwirkung 
zugesagt, während erfreulicherweise stetig neue Anerbieten zur Mitarbeit 
an der Sammlung einlaufen, worüber in meinen ,^tteilimgen'^ fortlaufend 
berichtet werden wird (die bereits erschienenen Bände sind mit zwei **, 
die unter der Presse befindlichen mit einem * bezeichnet): \ 

'^P. Baohmann, niedere Zahlentheorie. (Band X der Sammlung.) 
' M. Böoher, über die reellen Lösungen der gewöhnlichen linearen 

Differentialgleichimgen zweiter Ordnimg. 
. G. Bohlmann, Yersichemngsmathematik. 

G. n. Bryan, Lehrbuch der Thermodynamik. 

G. Castelnuovo, imd F. Enriques, Theorie der algebraischen Mächen. 
**B. Cznber, Wahrscheinlichkeitsrechnung und ihre Anwendung auf Pehler- 

ausgleichung, Statistik und Lebensversicherung. (Band IX.) 
*^Ij. E. Dickson, Linear Groups with an exposition of the Galois Field 

theorj. [In englisoher Sprache.] (Band VI.) 

F. Dingeldey, Kegelschnitte und Kegelschnittsysteme. 

F. Dingeldey, Sammlung von Aufgaben zur Anwendung der Differential- 

und, Integralrechnung. 

G. Eneström (in YerbinduDg mit andern Gelehrten), Handbuch der 

Geschichte der Mathematik. 

F. Enriques, Prinzipien der Geometrie. 

Flu Furtwängler, die Mechanik der einfachsten physikalischen Apparate 

imd Versuchsanordnungen. 
**A. Gleichen, Optische Abbildungslehre u, Theorie der optischen 

Instrumente* (Band VIII.) / 

M. Grübler, Lehrbuch der hydraulischen Motoren. 
J. Harkness, elliptische Funktionen. 
L. Henneberg, Lehrbuch der graphischeA Statik. 
SL Heun, die kinetischen Probleme- der modernen Maschinenlehre. 

G. Jung, Geometrie der Massen. 
G. Kohn, rationale Kurven. 

**A. Exaser, Handbuch der Lehre von den Thetafunktionen. (Band YII.) 
H. Lamb, Akustik. -AU- 

!B. V. Idlienthal, Differentialgeometrie. 



A, laoewy, Vorlesungen über die Theorie der linearen.Substitutipnsgruppen. 
**Q. laoria, spezielle, algebraische und transcendente Kurven der Ebene. 
Theorie und Geschichte. (Band V.) 
A. E. H. laove, Lehrbuch der Hydrodynamik. 

A. E. n. laove, Lehrbuch der Elasticitäi 

B. Mehmke, überi graphisches Rechnen und über Bechenmaschinen, 
sowie über nmneHsches Bechnen. 

W, Meyerhofer, die mathematischen Grundlagen der Chemie. 
**B, Netto, Lehrbuch der Combinatorik. (Band VII.) 
W. F. Osgood, allgemeine Funktionentheorie. 

E. Ovasza, aus dem Gebiete der Mechanik. 

**B, Pascal, Determinanten. Theorie und Anwendungen. (Band III.) 
S.' Finoherle, Funktional -Gleichungen und -Operationen. 
Fr. Fookels, Erystalloptik. 

A. PringsheinL, Vorlesungen über Zahlen- und Funktionenlehre. (Ele- 

mentare Theorie der unendlichen Algorithmen und der analytischen 
Funktionen einer komplexen Veränderliches.) Bd. L Zahlenlehre. 
Bd. n, Funktionenlehre. (B^nd I.) 
O. ^egre, Vorlesimgen über algebraische Geometrie, mit besonderer 
- Berücksichtigung der mehrdimensionalen Bäume. 

D. Seliwanoff, Differenzenrechnung. 

M. Simon, Elementargeometrie. ^ 

F. Stäröl^el, Lehrbuch der allgemeinen Dynamik. 

P. Stäokel, Differentialgeometrie höherer Mannigfaltigkeiten. 
O. Staude, Flächen und Flächensysteme zweiter Ordnimg. 
**0. Stolz und J. A. Gmeiner, theoretische Arithmetik, (Band IV.) 

B. Sturm, Theorie der geometrischen T^erwandtschaftenj 
!B. Sturm, die kubische Baumkurve. 

H. E. Timerding, Theorie der Streckensysteme imd Schrauben. 

E. Th. Vahlen, Geschichte des Fundamentalsatzes der Algebra. 
SL Th. VaMen, Geschichte des Sturmschen Satzes. 

A. Voss, Prinzipien der rationellen Mechanik. 
A. Voss, Abbildung und AbwicHimg der krummen Mächen. 
J. iG. 'Wällentin, Lehrbuch der theoretischen Elektrizität. 
**B. V. Weberi Vorlesungen über das Pfaffsche Problem und die Theorie 
der partiellen Differentialgleichungen 1. Ordnung. (Band 11.) 
*A, G. "Webster, the Dynamics of Particles, of rigid, elastic. and fluid Bodies 
beingLecturesonMathematicalPhysics. [inengUsoherSpUche] (Band XI.) 
A. Wiman, endliche Gruppen linearer Transformationen. 
W, Wirtinger, al)gebraische Funktionen imd ihre Integrale. 
"W. Wirtinger, partielle Differentialgleichimgen. 
n. G. Zeutben, die abzählenden Methoden der Geometrie. 

SS* Mitteilungen über weitere Bände werden baldigst folgen. 
Leipzig, Poststra&e 3. H (i TaiiVitiat* 

Oktober 1902. •*?• ^*« XÖUUllöl • 
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Vorwort. 



Die älteren Werke über Differenzenrechnung besitzen nicht die 
Strenge, die unsrer Zeit eigentümlich ist. In dieser Beziehung ist 
befriedigend das Buch von A. A. Mark off (Differenzenrechnung, 
St. Petersburg 1889—1891; deutsch von Th. Friesendorf und 
E. Prümm. Leipzig, bei B. G. Teubner. 1896). Ich habe versucht 
einige Verbesserungen zu machen. Vielleicht ist es mir gelungen 
einige Beweise zu vereinfachen und den Stoff systematischer an- 
zuordnen. Um die Arbeit den Studierenden zu erleichtem, habe ich 
speziellere Probleme beiseite gelassen und mich hauptsächlich auf 
elementare Fragen beschränkt. 

Worin nach meiner Meinung das Wesentliche in der Differenzen- 
rechnung besteht, habe ich in einem Referat dargestellt. (Encyklo- 
pädie der mathematischen Wissenschaften, Bd. I, p.918 — 937.) Dort 
findet man hinreichende Literaturangaben. In kleinem Räume konnte 
ich nur Sätze aussprechen, die Beweise sollte der Leser selbst suchen. 
Um ihm die erforderliche Mühe und Zeit zu ersparen, folge ich 
dem Vorschlag des Herrn B. G. Teubner und übergebe der Öffentlich- 
keit dieses Lehrbuch. 

St. Petersburg, den 6. Dezember 1903. 



D. Seliwanoff. 
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Erster Teil. 

Differenzen. 



Erstes Kapitel. 

Hauptsätze über Differenzen. 

1. Definitionen. Wenn eine Funktion f{cc) gegeben ist^ so 
heißt der Ausdruck f(x + h) — f{x) Differenz von f{x), in Zeichen: 

fix + h)- fix) = Afix). 

Wenn man diesen Ausdruck in eine andre Form bringt, so 
entsteht die Gleichung 

fix + h)-fix) = ipix) . 

oder 

Afix)^g>ix). 

Hier kann x beliebige Werte annehmen, h soll aber konstant 
bleiben. 

Ist z. B. fix) = x^y & = 1, so hat man 

^m = 1, A/-(2) = 7, Af{S) = 19, A/-(4) = 37, . . . 

Die Differenz von ^f{x) heißt zweite Differenz von f(x) oder 
Differenz zweiter Ordnung: 

AY(a:) = A/(a; + Ä) - A/"(«). 
Entsprechend ist 
A<^f{x) = AY(a; + Ä) - L^f(x), AYC«) = AYC« + ä) - AY(a;), . . . 
In dem Beispiel ist 

AY(0) = 6, AV(1) = 12, AY(2) = 18, . . . 
AY(0) = 6, AY(1) = 6,... 

AY(0) = 0,... 

Seliwanoff, Differenzenreolinnng. 1 
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Diese Resultate lassen sich in folgender Tabelle zusammenstellen: 
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2. DUterenzen einfooher Funktionen. Unter einfachen 
Funktionen versteht man diejenigen^ deren Differenz direkt bestimmt 
werden kann. Die Funktion heißt zusammengesetzt^ wenn ihre 
Differenz durch die Differenzen andrer Funktionen ausgedrückt 
werden kann. 

In der Differentialrechnung ist 

Die Differenz von a^ hat keine einfache Form, weil 



1 • A 



H-Ä" 



ist. Die einfache Funktion, die der Potenz entspricht, ist in der 
Differenzenrechnung die folgende 



f{x) = x{x — li){x — 2Ä) . . . (rr •— w — 1&). 
Die direkte Bechnung ergibt 



f{x + Ä) — f{x) =^(x + h)x(x — h)...(x — n — 2h) 

— x(x -'h)(x — 2h) . , .(x — n—lh). 

Um diesen Ausdruck zu transformieren, setzen wir vor die 
Klammem den gemeinsamen Faktor 

x(x — h)(x — 2h) . . .{x — n — 2h) 

und machen Reduktionen in den Klammem. 
Auf diese Weise findet man die Formel 



A[x(x — h)(x — 2h) . . .(x — n — lÄ)] 



= nhx(x — h)(x — 2h) . . . (a? — w — 2Ä). 

Spezielle Fälle dieser Formel sind 

Aa; = Ä, ^[x{x — h)] = 2hx, 

A[x(x — h)(x — 2h)'] = 3hx(x — Ä), . . . 
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Etwas einfacher ist die Formel 



. x(x — Ä)(aj— 2h) , , .(x — n—lh) , x{x — h){x — 2K) . . . (a? — n — 2ä) 

^ l-2.S...n "" 1.2.3...(n-l) ' 

die ganz analog bewiesen wird. 

Wir werden noch folgende einfache Funktionen angeben: 

9 m*, sinic, cosrr. 



x{x-^h){x-\-2h). . .(ic + n— lÄ) 

Hier ist m eine Konstante. 

Die direkte Rechnung ergibt 

1 



aj(a; + Ä)(a? + 2Ä)...(a; + n-lÄ) (a? + Ä)(a? + 2Ä)(ic + 8A) . . . (aj + nÄ) 

1 



a;(a: + Ä) (a; + 2Ä) . . . (a; + n- lÄ) 

L r-i--ii 

(a; + Ä)(a? + 2Ä) . . . (a; + n-lÄ)L« + «Ä a;J 

= — WÄ • aj(a. +Ä)(a? + 2Ä) . . . (x-^nh) 
Spezielle Falle dieser Formel sind 



x x{x-\-h) x{x-\-h) x{x-^h){x-\-2h) 

Etwas symmetrischer ist die Formel 

. l'2-8...(n— 1) , 1'2»3. ..n 

x{x-^h){x-^2h) . . . (x + iT^lh) "^ * x(x + h) {x + 2h) . . . {x + nh)' 

Eine sehr einfache Funktion ist m'^ weil 

Am* = in*+* — m* = m^(m^ — 1). 

Um die Differenz dieser Funktion zu bilden, genügt es, m* mit 
dem konstanten Faktor m* — 1 zu multiplizieren. 

Nach den bekannten Formebi der Trigonometrie ist 

Bixi(x + h) — siax = 2 sin — cos (a; + ö^) = 2sin — sin fa? + «^ + ^)^ 

cos(a: + h) — cosx — — 2 sin — sin (^•+-0) = 2sin — cosfrr + — + ^y 
Also wird 

Asina:= 2sin— cosfa: + -0] = 2sin — sin (rr H — 5—)' 

Acos.a? = — 2sin — sin (x + —j = 2sin ^^cosfa: + ^ )' 



1 
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Die wiederholte Anwendung der angegebenen Formeln liefert 
uns folgende Ausdrücke für Differenzen höherer Ordnung 



^ 1.2. 3...« '^ 1.2.8...(n-m) «»''»<♦»; 

y. n x(x — h)(x^2h),..(x — n-'lh) , n 

^ 1.2.3...n "* ^ 

. m l-2.8...(n-l) ^ , ,vm 1 • 2 ■ 8 . . . (n + m~ 1) 

a;(a; + Ä)(a?+2Ä)...(a; + n-lÄ) ^ "^ a;(Ä + Ä) . . . (a; + n + m-lÄ)' 

A~m'=m*(m*-l)-, 

A**sina; == (2 sin— j • sin (a; + w • ~2 )' 

A"cosir = (2 sin— j • cos (a? + w • ^ l' 

Das sind die wichtigsten Formeln^ die für das Folgende aus- 
reichen. 

3. Differenzen zusammengesetzter Funktionen. Wenn 

wobei C eine Konstante ist; so ist 

f{x + Ä) - f(x) = [g>{x + h) + C]- [ip(x) + C]=.q,(x + h)- q>(xl 

also £^f(x) = A(p(x). 

Bei der Berechnung der Differenz kann also der konstante Sum- 
mand weggelassen werden. 

ist, so ist 

£^f(x) = Cq>(x + h)- Cq>{x) = Clq>(x + Ä) - q>{x)'\ 

oder A [Cfy («)] = A ?> (x). 

Der konstaate Faktor kann also vor das Zeichen A gesetzt 
werden. 

Ist die gegebene Funktion f{x') Summe von mehreren Funktionen: 

f(x) = f,(x) + f,{x) + ••' + Ux), 
so findet man 

£^m=\f,{x+h)+f,{x+h) + -- .+f^(x+h)]-[f,(x)+f,(x)+-- ■+Uxy] 

=[A(«+Ä)-/i(a;)] + lf,(x + h)- f.ix)-] + ...+[f^(x + h) -f^ix)]. 

Es besteht also die Formel 

^[fi(^) + f%(-») + • • • + U{x)1 = A/i(a;) + A^,(a;) + • • • + LU{x). 
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Die bewiesenen Sätze lassen sich in einer Formel yereinigen. 

wobei C, C^, C^, , . ., Gfi konstant sind. 

4. Differenzen ganzer Funktionen. Es sei 

Es ist bequemer^ diese Funktion zuerst in die Form 
/ W = -^0 -r Al— 1 r^ 1T2 '"^8 TTTi ' 



. («--«)(«; — a — Ä).. .(ic — a — n — lÄ) 
•••+^ 1.2.3...n 

ZU bringen. Zur Bestimmung der Koeffizienten bemerken wir^ daß 
bei der Division von f{co) durch rr — a der Rest A^ ist und der 
Quotient 

Ai W = A + A 1.2 + ^ TTTs +' * ' 



, - (ic — a — Ä) . . .(o; — a~n— lÄ) 
••• + ^ 1.2.3...n 

Es ist noch zu bemerken^ daß f^ix) die Form hat 

/i(^)=l>o^"^+--- 

Dividiert man fiix) durch a; — a — &, so findet man den Rest A^ 
und den Quotient 



_r / ^ A 1 , A x—a — 2h , , M (a; — a — 2Ä) . . . (rc — a — w — lÄ) 
^8(a')=Al72+4» 1.2.3 +- + ^- l.2-3...n 

Das erste Glied von /^(a?) ist jPoüf*""*. 
Schließlich ist 



/• / \ ii 1 • ▲ x — a—n — lh , 

/^-»(^)=^-^- 1.2...(n-l) +^ 1.2-3...n i'oa' + ä- 

Hier ist g eine Zahl, die sich vollständig bestimmen laßt. 

Durch Division von /«— i(rr) durch x — a — n — 1 findet man 

den Rest ^ ^ **. — tt und den Quotient ^r-^ Dieser Quotient ist 

1 • 2 . . . ^w — 1^ 1 • a . . . n 

aber ^q. Es ist also 

Nachdem die Koeffizienten -4^, -4^, -4^, ...,-4n bestimmt sind, 
erhalt man mit Anwendung der Formeln (Nr. 2): 
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Af(x) ^hlÄ, + A^^-+''' + An i.2.3...(n-l) ^J' 

AY(a:X=Ä«[^ + ^,-^+... + Ä. 1.2.3.. .(n^2) ^> 

oder 

AY(ir) = r.l.2.3...w.i>o. 

Es besteht also der Satz: 

Die erste Differenz einer ganzen Funktion n^^ Grades 
ist vom (w — 1)**^ Grade und die n^ Differenz ist konstant. 

Für x = a ist 

Wir seheU; daß die KoefGlzienten A^, A^, . . ., An sich durch 
Differenzen ausdrücken lassen: 

Diese Werte der Differenzen bestimmen sich durch die Werte 
der Funktion f(x): 

f{a), f(a + h), f(a + 2h),...,f{a + nh). 

Da aber f(x) eine beliebige ganze Funktion w*®^ Grades ist, so 
besteht die Formel 



"^ 1-2-3.. .n~ * ~~^« 

für jede Funktion von der Form 

f(x)=pQaf' +PiX''-^ + p^af-^-\ \-pn-iX + pn, 

wobei einige Koeffizienten gleich Null sein können. 
Das ist die Newtonsche InterpolationsformeL 
Es sei z. B. 

a = 0, Ä = 1, f(x) = a;*. 

Wir haben schon in Nr. 1 gesehen, daß 

f(P) = 0, A/(0) = 1, AY(0) = 6, AY(0) = 6. 
FolgUch ist 

x^=x + dx(x — 1) + ie(x — l)(x — 2). 
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Als zweites Beispiel nehmen wir 

a==0, Ä = l, f(x) = l.2.3...n 

Hier ist y konstant vorausgesetzt. 

Nach der bewiesenen Formel (Nr. 2) ist 

^ /W— 1.2.3...(n-w) 

Daraus folgt, daß 

(a; + y)(a?4-y~l)...(a;4-y-n4-l) ^ y(y-l)...(y-n+l) 
1.2- 8 . ..n 1.2-3 .. .n 

j^x_ y(y-l)...(y-n + 2) a;(a?--l) y(y- 1) . . . (y-n + 8) 
' l' l-2.3...(n-l) "^ 1-2 * 1-2. ..(n — 2) "^ * * 



+ 



(n-1) * 1-2 1-2. ..(n — 2) 

x{x — l) . . . (rc — n-f 1) 



l-2-3...n 

Diese Relation besteht für beliebige Werte von x und y; sie ist von 
Cauchy benutzt worden, um den binomischen Lehrsatz auf beliebige 
reelle Exponenten auszudehnen.* 

6. Differenz des Produktes zweier Funktionen. Wenn 
man die Differenzen zweier Funktionen q){x) und i){po) kennt, so ist 
es leicht, die Differenz ihres Produktes zu bilden. . . • . 

Wenn f{(c) = (p(x)' ilf(x) ist, so ist 

Af(x) = q)(x + h)* tl}(x + h) — (p(x)' t(x) 

= [qp(^) + ^(p(poy]t{oc + h) - q)(x) • ilf{x) 
= v(o(^)[t(oo + h) — tipo)] + f(x + h) • Aqp(a?). 
Man findet also die Formel 

Z.B. für Ä = l ist 

6. Differenzen von sin u und cos u. Wir wollen voraussetzen, 
daß u eine Funktion von x ist, der^n Differenz bekannt ist. Es ist 
Asint« und Agosu zu berechnen. 

Wenn man x durch x + h ersetzt, so geht u ia u + Au über. 
Deswegen ist 

Asmw = sm(w + Au) — smM = 2sm— • cos(w +"0")' 
AcosM = cos(w + Am) — cosm = — 2sin-2^ • sin (u +-2~)* 
* Cauchy, Cours d'analyse, p. 159. Paris. 1821. 
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Damit beschließen wir die Sätze über Differenzen zusammen- 
gesetzter Funktionen. 

7. Berechniing Aquldlstanter Werte einer ganzen Funk- 
tion. Es sei f{x) eine ganze Funktion n^^ Grades^ deren Werte 

/'W; /*(ö^ + Ä), /-(a + 2Ä), . . ., f{a + nh) 

ausgerechnet seien. 

Es sollen die Werte 



/•(a + w + lÄ), f{a + n + 2h), /•(« + w + 3ä), . . . 
und 

f(a-h), f(a-2h), f{a-Sh),.,. 

bestimmt werden. 

Diese Aufgabe laßt sich leicht mit Hilfe der Differenzen auf- 
lösen. 

Durch wiederholte Subtraktionen bestimmt man die Reihe der 
Werte 



A f(a), A f(a + h), A f(a + 2h), . . ,, Af(a + n-lh), 
A* f{a), A» f(a + Ä), AY(a + 2h), . . ., AY(a + iT^Ä), 



n— 1/./ \ An— 1 



A-Y(a), A"-Y(a + Ä), 

Da aber l^f(x) konstant ist (Nr. 4); so ist 

ATfia -h Ä) = AY(a). 
Mit Benutzung der Formel 

A*-V(^ + Ä) = A*-V(^) + AV(^) 
und durch wiederholte Additionen erhalt man die Werte 

AT-'^fia + 2h) = A**"" V(a + Ä) -f A** f(a + h), 

A"""V(<^ + 3Ä) - A'^-VC» + 2Ä) + AT^^fia + 2h), 



A* /•(a-fn-lÄ)==A» /•(« + w - 2Ä) + A» /•(a-fw-2Ä), 

A f(a + nh) =A /'(a-f w-lÄ) -J- A* /'(a-f n~lÄ), 

f{a + n+ lÄ) = /*(a-fwÄ) +A /"(a-fwÄ). 

In derselben Weise werden auch die Werte 



f(a + n + 2h), f(a + n + 3h), . . . 
berechnet. 
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um f{a — h) zu bestimmen, bemerken wir, daß 

A7(o - Ä) = A"/-(o) 
ist, und daß die B«Iation 

A*-V(a;) = A*-V(a; + Ä) - A^f(x) 
besteht. 

Wiederholte Subtraktionen geben nns die Werte 

A"- V(a - Ä) = A"~ V(«) - ^" /■(« - Ä), 
A"- V(a - Ä) = A— V(«) - A"- V(o - Ä), 

A* f{a - Ä) = A« /•(«) - A» /•(» - Ä), 

A /"(a - Ä) = A /•(») - A* /•(» - h), 

f(a -h)= f(a) - A f(a- h). 

Ebenso findet man aach die Werte 

f(a-2h), /(a-3Ä), ... 
Als Beispiel nehmen wir 

a = — 1, Ä = l, f(x) = x^+x^-2x-l. 
Die B«chnm^ gibt 

fi-i)=i, m — 1, ^(1) — 1. 

Durch Subtraktionen findet man 

A/^C-l) 2, A/-(0) = 0, 

AY(- 1) = 2. 
Nach der bewiesenen Formel (Nr. 4) ist aber 

AY(^) = 6 
für jedes x. 

Die Anwendung der Additionen gibt uns 

AY(0) = 8, A/'Cl) = 8, f(2) = 7; 
A*/(l) = 14, A/(2) = 22, /(3) = 29. 

Um /"(— 2) und /"(— 3) zu bestimmen, sind folgende Rechnungen 
auszuführen: 

AY(-2) 4, A/-(-2) = 2, /•(-2) = -l; 

AY(-3) 10, A/'(-3) = 12, /■(-3)^ 13. 
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Diese Resultate lassen sich in der Tabelle zusammenstellen: 



X 


/ 


A/- 


AV 


AY 


3 


13 


12 


-10 


6 


2 


— \ 


2 


-4 


6 


1 




2 


2 


6 










8 


6 


1 




8 


14 




2 




22 






3 


29 






( 



Aus den Vorzeichen von f(x) erkennt man, daß die Gleichung 
f(ßi) = drei reelle Wurzeln hat, die in den Intervallen (— 2, — 1), 
(-1, 0) und (1, 2) liegen. 

Diese Methode fahrt immer zur Absonderung der Wurzeln einer 
numerischen Gleichung, sobald der absolute Betrag dßr Differenz 
zweier Wurzeln kleiner als Eins ist. 

8. Relationen zwischen den sukzessiven Werten nnd 
Differenzen einer Funktion. Wir stellen uns die Aufgabe, die 
Werte 

durch 

f{a), f(a + h), f(a + 2h), . . ., f(a + nh) 

auszudrücken. 

Um die Schreibweise abzukürzen, wollen wir folgende Bezeich- 
nungen einführen: 

f(a) = UQ, f(a + h) = Ui, /*(a + 2Ä) = w^, . . ., f(a + nh)=>Un^ 

Die gesuchte Relation kann auf symbolischem Wege ge- 
fanden werden. 

In der Gleichung 

Aujt = Uk+i — "Uk 

können symbolisch die Indizes durch Exponenten ersetzt werden. 
Wenn man aus der Gleichung 

den gemeinsamen Faktor u^ weglaßt, so wird 

A = w — 1. 
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Daraus folgt 

= u — mu H — - — U + ( — 1^ . 

Nach Multiplikation mit u^ und Ersetzen der Exponenten durch 
Indizes erhält man 

Diese Formel ist damit zwar gefonden^ aber noch nicht bewiesen. 
Das kann in verschiedener Weise geschehen. 
Wenn man in der Gleichung 

a durch a + h ersetzt, so vermehren sich alle Indizes um Eins und 
es wird 

Durch Subtraktion findet man 

In derselben Weise fährt man weiter fort 

A^Uk+i=Uk+s — 2uk^2+ w*+i, 

A^Uk = (Uk+Z — 2Uk+2 + Uk^i) — (Ma+2 — 2w*+i + Uk) 
= Uk+B — 3 Wa:+2 + 5Uk+i — Ukj 

A*Wa + 1 = Wa4-4— 3%+8+ 3wa4-2— w*+i, 

A*W* =(w*+4— 3Wa+8+ÄMa+2— W^+i) — (Wa+8 — 3W;fe + 2 + 3M*4.1 — Wjt) 
= W*+4 — 4wa+8 + 6^*4.2 — 4wa+i + Wifc. 

Damit ist die Richtigkeit der gefundenen Formel für 

m = 1; 2, 3 und 4 

bewiesen. Durch vollständige Induktion läßt sich zeigen^ daß die 
Formel allgemein gültig ist. 

Wir werden aber einen andern Weg einschlagen. Aus der 
früheren Betrachtung folgt, daß A^Uk die Form hat 

weil wir immer mit linearen und homogenen Funktionen der u zu 
tun haben. Die Koeffizienten A^, Ä^, . , ,, Am sind von der Funk- 
tion f{x)y von a und h unabhängig. Deswegen genügt es, diese 
Koeffizienten in einem speziellen Falle zu bestimmen. 
Wir setzen 

a = 0, Ä = 1, fix) = (f 
und findenr 



12 Erster Teil. Differenzen. 

Folglich ist 

(c— l)'« = ^(-»+ J^cm-i + ^(«»-a+ \.A^, 

Die Koeffizienten J. sind aber ans dem binomischen Lehrsatze 
bekannt. 

Es ist also allgemein 

A = l, A=-m, J,= ?!^^,..., ^.= (-1)- 

Setzt man in der bewiesenen Formel 

Ä; = 0, m = 1, 2, 3, . . ., w, 
so erMlt man die gewünschten Formeln 

A*Wo=Ws,— 2wi + Wo, 
A^Uq = Wg — 3wj, + 3i^ — Uq, 

An I ♦* (♦* — 1) , / -t\n 

UQ=Un—nUn^l+ '^ g ' Un-^i h (— 1) % 

8. Jetzt wollen wir znr umgekehrten Aufgabe übergehen. Es 
sollen die Werte 

f(a), fia + h), f(a + 2h),..,, f(a + nh) 
durch 

f(ä), Af{a), AY(a),..., AV(a) 

ausgedrückt werden. 
Aus der Formel 

folgt 

Symbolisch ist 

oder 



w*H-i = w*+ Am*. 
w*+^= M*+ Aw* 

M = l + A. 



Daraus ergibt sich 

w- =l + mA + ?^^?^A^+..:+A", 

w'»+*=w*+mAw*+?!^^^AV+... + AV, 

X * A 

I A i^(*'*~1)a« I lAin 

«^m+*=w*+ wAwi+--Y-2-^A*M*H hA Uk. 
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Um diese Formel zu beweisen, beginnen wir mit den Fällen 
m = l, 2, 3. 

Wir haben schon gesehen, daß 

Daraus folgt (Nr. 3) 

Auk+i=£^Uk+ A^Uk, 

Uk+2={uk+ Auk) + (Awa+ A*Wa)==Wa+ 2Awt+ A*Wi, 

Aw*+8= Awa+2A*w*+ A^w*; 

Uk+fi= (Uk+ 2Aw*+ A^Uk) + {AUk+2A^Uk+ A'w*) 

== wt+ 3Awt+ 3A*w*+ A^w*. 

Wir haben hier immer lineare und homogene Funktionen von 
Uk, Auk, A^Uk, . . .; folglich ist 

M*+„»= JB^Uk +BiAuk+ B^ A^Uk H h B^A^^Uk, 

wo die Koeffizienten von der Wahl der Funktion f(x) unabhängig sind. 
Für 

findet man 

6*+»»= 5oC* + B^(^{c - 1) + B^€^{c - 1)2 + ... + 5^c*(c - ly 
oder ^=5^+5^(^_i) + :B^(c_1)2+... + 5^(^_1)^^ 

Setzt man c = l + t ein, so wird 

(1 + ty=B^+ BJ + J5, <»+••• + £„<". 
Die gesuchten Koeffizienten haben also die Werte 

B, = l, B, = m, JB^ = ?^^L=i), . . ., :b, = i. 

Es bestehen daher die Formeln 
«*i = «^o+AWo> 



Un=-%+nAuQ+ ^^^ g^^ A* Wo + h A" Wo- 

Alle diese Formeln können durch die Formel 

umfaßt werden. Hier kann nur die Werte 

a«r,al.rr..r, ^ = 0, 1, 2, 3, . . ., W 

annenmen. 

Damit schließen wir dieses Kapitel und gehen zu Anwendungen 
der bewiesenen Sätze über. 
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Zweites Kapitel. 
Interpolation. 

10. Oenaue Interpolation. Die einfachste Aufgabe ist die 
folgende: Es sind die Werte einer ganzen Funktion f{() von nicht 
höherem als n^^ Grade gegeben 

f{a), f(a + h), f(a + 2h),..,, f(a + nh). 
Es soll der Wert der Funktion für t = x berechnet werden. Diese 

Aufgabe war schon gelöst fiir den Fall; daß , gleich einer ganzen 
Zahl ist (Nr. 7). 

Die bekannten Werte können wir wieder mit 

Uq, Ulf tig, . . ., Un 

bezeichnen. Es ist also der Wert der Funktion u, gegeben für 

jgf = 0, 1, 2, . . ., w. 

Nach der bewiesenen Formel (Nr. 9) ist 

I A I ^(^— l)Aa I I e(z—l)...(z — n + l)j.n 
W,= t*0+^At*o+-i-2-'A^Wo+--- + 1.2. 3.. .n ^ ^0 

für diese Werte von 0. 

Wenn wir setzen 

t—a 

SO ist 



+ 1.2.8...n ^ "o 

eine ganze Funktion, bei der die Relation 

für w + 1 verschiedene Werte von t besteht Diese Relation muß 
identisch erfüllt sein, weil die algebraische Gleichung F{{) — f{() = 
von nicht höherem als nf^^ Grade nicht mehr als n verschiedene 
Wurzeln besitzen kann. Wenn wir t = x setzen, so finden wir 

+ ^ 1.2.8. :.n -^ ^W- 

Das ist die Newtonsche Interpolationsformel, der wir schon 
früher begegneten (Nr. 4). 
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11. Andre Aufgaben der Interpolation wollen wir nur für ganze 
Funktionen von nicht höherem als dem 4*®^ Grade behandeln. 

Es sind Werte einer ganzen Funktion f(t) vom 4*®^ Ghrade ge- 
geben 

f{a), f{a + h), f(a + 2h), f{a + 3h), f(a + U). 

Man soll für t = x den Wert der Funktion und ihrer Ableitungen 
berechnen. 

Nach der Newtonschen Formel ist 

m = f(a) + ^ . Af(a) + ^^^t^V.^"'^ ^V(a) 

+ 1.2.3Ä» ^ fW + 1.2.8.4Ä* ^ / W 

für jeden Wert von t 

Der Einfachheit wegen setzen wir t = a + hz, die Punktion f(i} 
verwandelt sich dabei in F(0) und wir finden 

F(,) = f(a) + ^Af(a) + 4^ AY(«) + "^'L^V'^ ^'/"^ 



oder 



<i(g-l)(g-2)(g-8) . ^ , , . 
"^ 12.84, ^ ' W 



+ (^-6r'+ll.«-6.).5^. 
Durch Differentiation erhält man die Formeln 
F' (^) = A/'(a) + (2.-l).^ + (3.«-6. + 2).^ 

+ (4,«- 18.«+ 22. -6). 5^, 



JP" (.) = A'f(a) + (. - 1) AY(a) + (6«*- 18. + 11) • ^^, 

F'" (.) = AY(a) + (ßg - 3) ^^^, 

J""'(.) = AY(«). 

Wenn . = y für t = x ist, so sind die gesuchten Werte 

f'ix)=>iF'(y), f"(x)=j-,F"{y), f"{x) = ^r" iv), 

rix)^lF^(y). 



x — a 



Hier ist y durch —j— zu ersetzen. 

Es seien z. B. für eine ganze Funktion f(t) von nicht höherenk 
als 3^ Grade die Werte 
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f(0)=l, A2)==3, r(4)==53, /•(6)-199 
gegeben. 

Es soll f"(l) berechnet werden. 
In diesem Falle ist 

/•"(l) = ^[AY(0)-i-AY(0)] 

AY(0) = 53 - 2 . 3 +1 = 48, 
AY(0) = 199 -3-53 + 3-3 -1 = 48, 

r(l) = ^(48-24) = 6. 

12. Für eine ganze Funktion 4^^ Grades nnd für A =» 1 seien 
die VT erte 

m, A/-(0), AY(0), AY(0), AY(0) 
gegeben. Es sollen die Werte 

berechnet werden. 

Wenn wir die Bezeichnungen 

dfia) = f(a + ^) - f(a), d'f(a) = df(a + ^) - 8f{a), . . . 

einführen^ so ist nach der Newtonschen Formel 

OT(w-l)(OT-2)(m-8) M ^ ,^v 
+ 1.2-3.4 *^^ W- 

Die Werte 8f((S), S^f(0), dY(0) und d*/"(0) lassen sieh ans der 
Formel 

m- m + 'Äfto) + '4t ''^■''(o) + "'7'^",""" 'i'fto) 

bestimmen, die für jedes t gültig ist. Man findet 
/■© = AO) + Ä Ar(0) - 4 AY(0) + i^ AY(0) - S ^V(0), ■ 

^(i^) = /"(0)+^^/^(0)-ä^V(0)+SAY(0)-S^VW, 

/■(Ä) = ^W + Ä ^^CO) - ä AY(0) + j^ AY(0) - ig^ AY(0). 
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Die Reihe von Subtraktionen gibt 

sm = ^ A/-(0) - 4 AY(0) + ^ AY(0) - ^ AY(0), 

*/•© = ro^m-iö ^V(0) + 3^ A»^(0) - S AY(0), 

^/^S) = ho^m - wo^'m + ä^^V(0) - 8^ AY(0). 

*V(0) = i^ AV (0) - jA_ A»;.(0) + ^ AY(0), 
*V© = iSö AY(0) - ^ AY(0) + ^ AY(0), 
*V© = i^ AY(0) - ^ A»AO) + ^ AY(0). 

*'/'(0) = i^AY(0)-^AY(0), 

*vQ = iÄö^V(0)-2^AY(0). 

Die gesuchten Formeln sind also 

*/(0) = ^ A/^(0) - ^ AY(0) + ä^AY(0) - S^V(0), 

*V(0) = iJö AY(0) - j^ AY(0) + ^ AY(0), 
*V(0) = i^AY(0)-^AY(0), 

Wir haben schon gesehen (Nr. 7), daß die Gleichung 

im Interyall (-- 2, — 1) eine Wurzel besitzt. Die Gleichung 

f{x) = 9 (- 2 + oj) = 
hat aber eine Wurzel zwischen und 1. 



Seliwanoff, Differensenreohnong. 
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Nach der angegebenen Methode können die Werte 

berechnet werden. 

Wir wollen voraussetzen, daß /*(^) tind f[ ^^ ) verschiedene 
Vorzeichen haben. Dann hat die öleichnng (p(t) = eine Wurzel 
im Intervall (- 2 + ^, - 2 + ^). 

Jetzt ist in q){t) die Substitution 

t 2 + ^ 

anzuwenden. Wir erhalten eine neue Funktion 



/iW = 9'(-2 + *-^) 



die auch eine Wurzel zwischen und 1 besitzt. 

Nach demselben Verfahren findet man zwei Zahlen ~ und ' ^ > 
für die /ife) und A T \q ) verschiedene Vorzeichen haben. Die 
Wurzel der Gleichung qp (<) = liegt also im Intervall 

(-^ + ro + w' -^ + ro + '4^y 

Auf diese Weise kann man die Wurzel der Gleichung mit beliebiger 
Genauigkeit berechnen. 

13. Wir wollen jetzt folgende Aufgabe der Interpolation losen. 
Es ist eine ganze Funktion 4*®^ Gh^des zu bilden, die für 

/•(«), /"(«), /■(« + Ä), r{a + h), f(a + 2h) 

vorgeschriebene Werte besitzt. 
Die Funktion 

Fi^) = m + ^A^(«) + (^-«)(^^-«-^) AV(a) 

hat die Eigenschaften 

F(a) = /•(«), F{a + h) = f{a + h), F(a + 2h) = f(a + 2h). 

Dieselben Eigenschaften bestehen atacb. fOr 

*(«) = F{x) + (x — a)(x — a — h)(pi; — a — 2h)tp(x), 

wie aach die Funktion qj(x) gewählt wird. 
Aus den Fordenu^en 

*'(a) = f'(a), 9'(a + Ä) = f'(a + h) 
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lassen sich die Werte (p{a) und ^(a + Ä) berechnen. Deswegen ist 
die Funktion 

Yollst'ändig bestimmt. 

Es sollen z. B. die Bedingungen 

/(o)=i, /'(o)=i, m)=2, /"(i)=3, m=^i 

befriedigt werden. In diesem Falle ist 

A/'(0) = 1, A/'(l) = 5, AY(0) = 4, 

F{x) = 1 + a; + ^^^.4 = 2a;»- a; + 1, 

«(a;) = 2a:* - a; + 1 + a;(a; - 1) (a; - 2)(p(x) 
= 2a;» - a; + 1 + (a;» - 3«* + 2a;)9>(a;), 
*'(a;) = 4a; - 1 + (a;» - 3a;« + 2a;)9'(a:) + (3a!* - 6a; + 2)y (a;), 
«'(0) = - 1 + 2<p(0) = 1, 9(0) = 1, 
*'(l) = 3-y(l) = 3, 9(1) = 0, 

9(ir) «= 1 + a; • (— 1) = — i» + 1; 
f{x) = 2x^-x + 1 + (x^'-Sx^+ 2x)(- x+1) 

^-a^+4x^-3x^+x+ 1. 

Wir wollen noch ein Beispiel behandeln. Es soll eine ganze 
Funktion 3*®"^ Ghrades die Bedingungen erfüllen 

/•(0) = 1, ^(1) = 1, f{l) = l, /"'(1) = 4. 

Nach dem Taylor sehen Satze hat die Funktion 

F{x) == 1 + (äj - 1) + 4^^^ =^ 2x'^ 3a; + 2 

die Eigenschaften 

JP(1)-1, JP'(1) = 1, i^"(l) = 4. 

Dasselbe gilt auch für 

0{x) = 2x^- 3a? + 2 + (a; - 1)V(^) 

bei willkürlicher Wahl von q>(x). Da aber nur eine Bedingung zu 
erfüllen ist^ so können wir q>(x) konstant voraussetzen. Es soll 
©(0) = 1 sein, folglich ist 

2-«)p(0) = l, y(0)=l 

^^ /•(a;) = 2a;*-3a; + 2 + (a;-l)»=a;»-a;*+l. ., 

Dieselbe Aufgabe konnte man in andrer Weise lösen. 

2* 



* j - ■> j 
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Den Bedingungen f(P) = 1, f(l) = 1 genügt die Funktion 

0(x) = 1 + x(x — l)9(i») = 1 + (a;* — x)q>(x). 
Die Ableitungen dieser Funktion 

0f(x) = {x^-'X)g)' (x) + {2x - l)g){xl 
Q'\x)=^(x^- x)q>"(x) + 2(2a; - l)(p'(x) + 2(p{x) 
haben für x = l die Werte 

0\1) = (p(l) = 1, ^''(1) == 2q>\l) + 2q>(l) = 4. 

Daraus folgt 

9)(1) = 1, 9>'(1) = 1, 

q)(x) = 1 + (a? — 1) = a?, 

f(x) = 1 + (x^— x)x = a?* — a?* + 1. 

Granz ähnlich kann man ganze Funktionen höheren Grades bilden^ 
die vorgeschriebene Werte 

/(a + Ä), f{a-V Ä), r (a + Ä), . . . , /"K-D (a + Ä), 

/•(a + mh)y f{a + mÄ), /•"(« + mÄ), . . ., f^''m-^){a + wä) 
besitzen. Der örad dieser Funktion kann nicht die Zahl 

a+ai+Os + '-'+Om— 1 

übersteigen. (Vgl. Markoff, Differenzenrechnung, S. 1 — 6.) 

In den betrachteten Fällen der Interpolation wurde vorausgesetzt, 
daß die Werte der unabhängigen Yariabeln in äquidistanten Intervallen 
auftreten, weil nur in diesem Falle die Differenzen anzuwenden sind. 

14. Angenäherte Interpolation. Für Werte von i, die in 
einem Intervall liegen, wollen wir die gegebene Funktion f(t) an- 
genähert durch die ganze Funktion F{t) ausdrücken, die gewissen 
Forderungen genügt. Der Fehler dieser Annäherung soll abgeschätzt 
werden. 

Z. B. soll die ganze Funktion F(t) vom nicht höheren als 
^ten (Ji-ade sein und die Gleichung 

erfüllen für 

^ = a, a + Ä, a + 2Ä, ...,a + 



• • Es ist uns schon bekannt, wie man eine solche ganze Funktion 
i'Ät (Nr. XO). 
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Es sei X eine ZaU; die im gegebenen Intervalle liegt. Wir 

werden nähemngsweise 

f{x) = Fix) 

setzen und es ist der Fehler dieser Annäherung abzuschätzen. 
Die Punktion 

fit) - F{t) -- K i.2.3...(n+l) 

bei beliebigem Werte von K verschwindet für 

t^a, a + h, a + 2h, . . . , a + nh. 

Wir wollen K so wählen, daß diese Punktion noch für t = x 
verschwindet. Hier ist x eine Zahl, die von a, a -{- h, a + 2ä, . . . , 
a + nh verschieden ist. 

Die Gleichung 

*(0 = m -nO - ^ i.2.3...(n + l) ^ 

hat dann mindestens n + 2 Wurzeln. Nach dem Roll eschen Satz 
hat die (n + 1)*® Ableitung von ^(t) mindestens eine Wurzel, die 
wir mit | bezeichnen wollen. Diese Zahl ist eine mittiere zwischen 
der größten und der kleinsten unter den Zahlen a, a + nh und x. 
Es ist aber 

folglich ist 

Es besteht also die Pormel 



.(x—a){x—a—h). . .(x—a—n—lh) ^nn/ \ , (x-a)(x-a-h), . ,(x-a- nh) r^n^t)rt\ 
+ 1.2.8...n.Ä« ^ ^^*^+ 1.2.8...(n + l) f ^^j- 

Das ist die New ton sehe Interpolationsformel mit dem Rest- 
gliede von A. Cauchy (vgl. Cauchy, Oeuvres, 1 serie, tome 5, 
p. 422 oder Comptes rendus, tome 11 [1840], p. 787). 

Besonders sind die Pälle zu bemerken 

M = 0, f{x) = f{ä) + (x-a)f'{i), 
■ n = l, fix) = fia) + ^ Afia) + <.--)(---^'> ^''(|), . 

n = 2, fix) = fia) + ^« Afia) + ^"-f.^r.^"'^ ^V(a) 

{x-a)(x-a-h){x-'a--2h) fitift\ 
+ 17273 f W- 

Die erste dieser Pormeln ist die Lagrangesche. 
Durch Weglassung der Restglieder entstehen Näherungsformeln, 
deren Pehler durch die Restglieder ausgedrückt sind. 
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16. Wir wollen einen besonderen Fall der Interpolation betrachten. 
Für t = x soll die gegebene Funktion f({) angenähert durch eine 
solche ganze Funktion 3*®^ Gh»des F{t) dargestellt werden, daß die 
Bedingungen 

F(-l) = fi-l), F(0)=^fiO), F'(0) = f'(p), Fil) = f(l) 

erfüllt seien. 

Die Funktion F(t) hat die Form 

wo die Konstante Ä sich aus der Bedingung 

-F"(o)=r(o) 

bestimmen laßi Wir finden 

JF'CO = A/(- 1) + ^(2t + 1)AY(- 1) + {(3<»- 1)^, 

f(0) = Afi-l) + ^L'fi-l)-^Ä, 

A = QLf{- 1) + 3 AY(- 1) - 6/"(0), 

A = 6[/-(0) - A- 1)] + 31/(1) - 2m + /•(- 1)] - 6^(0). 

Nach Reduktionen ergibt sich 

A = 3/^(1) - 3A- 1) - 6^(0). 

um die Interpolationsformel mit dem Bestgliede zu finden, 
betrachten wir die Funktion 

q,(t)=m-F{f)-K'^^^^^^, 

wo die Konstante K so gewählt ist, daß 9{x) gleich Null ist. 
Wenn wir die vier Zahlen 

- 1, 0, 1, X 

der Größe nach ordnen, so entstehen solche vier Zahlen 

«1, ttg, aj, a^j 
die den Bedingungen genügen 

«1 < Og < «3 < a^. 

Diese Zahlen genügen der Gleichung ä^(^) = 0. Nach dem 
Roll eschen Satze hat die Gleichung 9^{t)^0 mindestens drei solche 
Wurzeln 6^, t^, tj, daß 

% < 6i < «2 < &2 "^ ^8 "^ ^8 "^ ^4- 



t ^ • • • 

= ~ • ♦ • 
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Keine der Zahlen b^y b^ und b^ ist also gleich Nnll. Außerdem 
ist diese Gleichung noch fiir ^ = befriedigt. Die Gleichung 0'(t)=^O 
hat also mindestens vier verschiedene Wurzeln. Nach dem Bolle sehen 
Satze hat die dritte Ableitung von 9\t) oder 9^{t) mindestens eine 
Wurzel < = I, die zwischen a^ und a^ liegt. Da aber 

Wir haben also die folgende Formel bewiesen 

f{x)=^f{^ 1) + ^ a;^(^ 1) + ^-gi^ AV(-^ 1) + ^ ^" v.^,"^r'^ 

Der Koeffizient A hat den angegebenen Wert. 

16. In andern Fällen werden die Bestglieder mit Hilfe der- 
selben Betrachtungen gebildet. 

Wenn z. B. die Näherungsformel 

genau ist für 

t = ay a-\-hy a -\- 2h, a + 3ä, 

f(f)^F\t) für t = a, a + h, a + 2h, 

f"(t) = F"({) für t=^a, 

dann ist F(t) eine ganze Funktion vom 7*®" Grade und die Näherungs- 
formel f(x) = F(x) enthält den Fehler 

(a: — a)»(rc — a — Ä)'(a; — a — 2Ä)»(a; — a — 3Ä) 



1.2-3-4-6-6-7-8 



f^'m^ 



Die Zahl I ist zwischen der größten und der kleinsten unter 
den Zahlen a, a + 3ä, x enthalten. (Vgl. Markoff, Differenzen- 
rechnung, S. 6 — 8.) 

17. Anweudnng der angeuftherteu Interpolation auf die 
Berechnnng der Wnrzeln nnmerischer CHeichnngen. Wir 

machen die Voraussetzung, daß die numerische Gleichung f(t) = 
im Intervall (a, 6) nur eine Wurzel besitzt und daß f"(t) immer 
dasselbe Vorzeichen in diesem Intervall behält. 

Da die Vorzeichen von f(a) und f(b) verschieden sind, so ist 

Wir werden f(t) naherungsweise durch eine solche lineare 
Funktion jF(^) ausdrücken, für die 

F(a) = f(d) und F{b) = f(b). 
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Hier ist h = b — a und deswegen ist 

F{t) = /•(«) + ^ 1/(6) - fia)-]. 
Für ^ = a? wird 

Den Wert rc wählen wir so, daß F(x) = wird. Polglich ist 

Es liegt also x zwischen a und &. Für diesen x kann nicht 
/*"(!) = sein. Die Gleichung 

wäre dann befriedigt für t = ayX und &. Nach dem Bolleschen 
Satze hätte die Gleichung ^(f) = mindestens zwei Wurzeln zwischen 
a und h, und die Funktion 

*"(0 = f"(t) 

müßte in diesem Intervall ihr Vorzeichen wechseln, was gegen die 
Voraussetzung ist. 

Weil F{x) = ist, so ist 

Der betrachtete Wert von x läßt sich einfach durch Deter- 
minanten ausdrücken. Es ist 



X = 



m f{b) 

a h 



m m 
1 1 



Wenn t im Intervall (a, 6) liegt und 

na)r{t)>o 

ist, so ist f{a)f"{l)>0 und folglich 

f{a)nx)<0, 

da {x — a){x — h) negativ ist. Die Wurzel der Gleichung f({) = 
liegt also zwischen a und x. 

Wenn wir a und b als erste Annäherungen der Wurzel be- 
trachten, so sind 

zweite Annäherungen. 
Es sei 

/•(a)r(0<0, oder fih)f"it)>0, 

für Werte von t, die im Intervall (a, h) liegen. 



-. -» — ^ * '^ - •H^aiBäHHff9>'^iOT" 
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Dann ist 

f{V)r{i)>0 und f(h)fix)<0. 

Die Wurzel von f{t) = liegt zwischen b und x. Die zweiten 
Annäherungen sind in diesem Falle 

In derselben Weise kann man von a^ und \ zu den folgenden 
Annäherungen Og und &2 übergehen usw. (Vgl. Weber, Lehrbuch 
der Algebra, I, § 116.) 

18. Anweudnng der angeuftherteu Interpolation anf die 
Berechnnng der Logarithmen nnd Antilogarithmen. Wir 

nehmen z. B. eine Tabelle fünfstelliger gewöhnlicher Logarithmen 
aller ganzen Zahlen von 1000 bis einschließlich 9999. Mit Hilfe 
der Werte log N und log {N + 1) soll log (-W+ oo) berechnet werden, 
wobei < a: < 1 ist. 

Die Literpolationsformel 

fix) = m +xLm + ^'V"(i) 

gibt die Näherungsformel 

log {N+ x) = log iV+ X [log {N+x)- log JV] 

mit dem Fehler 

x{l — x) log e 

Es ist 

6<3<>/iÖ, l0g6<J; 

Daraus folgt 

^<j^e oder 6<7.10-8. 

Dieser Fehler hat keinen Einfluß auf die fünfstelligen Logarithmen« 

19. Die Berechnung der Antilogarithmen führt auf die Frage: 
Es sind die Werte gegeben: 

log j!V'= a, log (iV+ 1) = l. 

Aus der Gleichung 

log(JV+y) = a + a; 

soll y berechnet werden für einen Wert von Xy der zwischen und 
6 — a liegt. 
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Wir setzen ^. . _,, . ^/v . 

Die Anwendung der Interpolationsformel 

gibt die Naherong 

mit dem negativen Fehler 

8 = _ ^(&-a-^) .iOa + f(l 10)8.* 

Es ist 

,(j_,_.)=(»^«)'_(-f_.)'^(?^)'. 

Ans der Relation 

folgt «^^" = 1^ 

110 = .— • 

log« 
Da aber 10'»+^<10> ist, so ist 

Wenn wir bemerken, daß 

10» = JV^+ 1, 6 - a = log (JV+ 1) - log JV= log(l + ^), 
so finden wir 

-a<^(N+ 1). (löge)-» [log (l +^)f 

Dieser Ausdruck kann noch, vereinfacht werden. 
Für positive z ist 

q)(/) = log (1 + is) — jsloge< 0, 

weU <p(p) = und 

sind. Deswegen ist 

log (l + ^) < i^ log ^ 
und folglich ^y+i 

Da aber N+l£10^, N^ 10», so ist die positive Zahl 

-6<2.10-8. 

In diesem Eapitel ist die Interpolation nur als Anwendung der 
Differenzenrechnung behandelt worden. Ausführliche Mitteilungen 
über diesen Gegenstand findet man im Referat von J. Bauschinger 
(EncykL der Math. Wiss., I, S. 799—820). 

* 1 10 ist der natürliche Logarithmus von 10. 
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Drittes Kapitel. 
Angenäherte Bereclinnng bestimmter Integrale. 

b 

20. Elnleiteude Bemerkungen, um das Integral 1 f(x)dx, 

a 

WO a < 6 ist, angenähert zu berechnen, gebraucht man eine der 
Interpolationsformeki (Nr. 14 — 16) 

Hier sind F(x) und q)(x) rationale ganze Funktionen und | ein 
Mittelwert. 

Nach Integration findet man die Näherungsformel 

b b 

Jf(x)dx=^ JF{x)dx 

a a 

mit dem Fehler 






a 



Es ist I eine unbekannte Funktion von X] wir wissen nur, daß 
für a<Cx<,b auch a < S < 6 ist. Der Fehler der Näherungsformel 
ist nur dann abzuschätzen, wenn die Funktion q)(x) für alle x zwischen 
a und 6 ihr Vorzeichen behält. Nach einem Satze der Integral- 
rechnung ist dann 



-r^Kv)fv(^)dx, 



wo a<rj <b ist. 

Die Aufgabe ist also auf die Integration ganzer rationaler 

Funktionen zurückgeführt. Die Rechnungen sind noch einfacher, 

wenn die Grrenzen des bestimmten Integrals — 1 und + 1 sind, z. B. 
1 






(j?o + Pi^+ -Pa^^ + Pü^ + '9i?^)äx = 2jPo + \'g^ + |i?4. 



Um die Grenzen des gegebenen Integrals in — 1 und + 1 zu 
verwandeln, setzen wir 

und es sollen die Bedingungen erfüllt werden 

^(-l) = a, ^(1) = 6. 
Nach der New ton sehen Interpolationsformel ist 

a? = a + -^(6 - a) = -^< + -^. 



28 Erster Teil. Differenzen. 

Yermöge dieser Substitution findet man 

b 1 

a —1 

In den folgenden Paragraphen werden wir verschiedene Inter- 
polationsformeln anwenden. Darans entstehen verschiedene Methoden 
zur Berechnung des gegebenen Integrals. 

21. Methode der Bechtecke. In der Formel 

b 1 



ffix)dx = ^-^ß{t)dt, 



a — 1 

wo 



setzen wir (Nr. 14) 

j^(t) = j^i- 1) + (t + 1) J^'(<i), (- 1< <i < 0- 

Daraus entsteht die Näherungsformel 



ß{t)dt = 2/(- 1) 



mit dem Fehler 

1 



Hier ist — 1 < r < 1. 

Um die Funktion f{x) einzufahren, bemerken wir, daß 

Wir finden die Naherungsformel 

b 



Jf{x)dx = (6 - a)fia) 



mit dem Fehler 

B = \{h-af.f\%), 

WO a < I < & ist. 

Geometrisch ist das gegebene Integral eine Flache, die in der 
Näherungsformel durch ein Rechteck ersetzt wird. Daher rührt der 
Name der Methode. 

Um den Fehler zu vermindern, zerlegen wir das Intervall (a, 6) 
in n gleiche Intervalle. Das Integral verwandelt sich in die Summe 

b a-\ a+2 b 

f{x)dx = / f{x)dx + 1 f(x)dx+''' + j f(x)dx. 

•y ^, b — a €/ -b—a 

a a a-\ a+« — 1 
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Kachdem wir auf jeden Smumaaden die Käherungsformel angewandt 

haben, finden wir 

b 



ß 



f{x)dx = ^ (jTo + ^1 + y» + • • • + Vn-l), 



a 

WO der Kürze wegen die Bezeichnung 

y* = /*(« + Ä • -^) 

eingeführt ist. 

Der Fehler dieser Näherungsformel ist 

«=v^'[/''(ii)+r(y + -"+r(i«)], 

wo a + (Ä — 1) -^^^ < I* < a + Ä; -^^^ ist. 

Diesen Ausdruck kann man vereinfachen. 

Es wird vorausgesetzt, daß f{x) endlich und stetig im Intervall 
(cby h) ist. Diese Funktion hat dann eine obere Ghrenze A und eine 
untere Grrenze J?.* 

Aus der Reihe von Ungleichheiten 

nA ^ /'(y + /•'(!,) + ... + /•'(!„) ^ nB 

Nach einem Satze von Weierstraß kann f{(c) jeden Wert 
zwischen A und B annehmen. Deswegen gibt es eine solche Zahl %, 

^"^^ r(s.)+r(i.)+---+r(i«) _ ^,/|)_ 

Es ist also 

Wenn wir die Interpolationsformel 

benutzen und dieselben Rechnungen wiederholen, so finden wir die 
Näherungsformel 



/' 



f(x)dx = ^(Vt + y, + • • • + jr„) 



a 

mit dem Fehler ^ ,, ,, 

* Tann er j, Introduction k la thäorie des fonctions d'une variable, 
p. 101 et 126. 



30 Erster Teil. Differenzen. 

Wenn f{po) immer dasselbe Vorzeichen im Interrall (a, &) hat, 



80 ist der Wert von 1 f{x)dx zwischen 



a 

b — a 



(vo + yi + yi + '-' + yn-i) 



n 
nnd 

-^(vi + ys + yz + '-'+yn) 

enthalten; weil s nnd s^ verschiedene Vorzeichen haben. 

22. um das Integral 

b 1 

Jf{x)dx = ^ß(t)dt 

a — 1 

zn berechnen^ kann man nicht die Formel 

gebrauchen^ weil t, der Paktor bei J^^t^), das Vorzeichen im Intervall 
(— 1, 1) wechselt nnd der Fehler der Nähemngsformel nicht ab- 
geschätzt werden kann. Es ist zu setzen (j,M.arkoS^, S. 53) 

und man findet 

1 



/ 



/(0d^ = 2J=(0) + {j="(r). 



3 

— 1 



DarauB folgt die Nähemngsformel 

b 

^a + 6> 



/' 



a 

mit dem Fehler 

24,(P ^) I VS;- 2 1.2.3 1.2 

um eine genauere Formel zu bekommen^ zerlegt man das Intervall 
(fly V) m n gleiche Intervalle. Wenn man dieselbe Methode wie in 
Nr. 21 benutzt; so findet man die Näherungsformel 

h 



/' 



f(x)dx = -^ (y£+ ys + h ya«-i) 



2 S 

a 

mit dem Fehler . ., .. .„,.. 

2n» * 1 . 2 . 3 * 1 . 2 * 

Hier ist wieder 

y,-f((^+^~} 
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28. Methode der Trapeze. Zur Bereclinimg des Integrals 

b 1 

ff(x)dx = ^fj'(t)dt 

a —1 

werden wir die Funktion S(t) angenähert durch, eine lineare Funktion 
^(t) ausdrücken^ die deil Bedingungen genügt 

«(-.1) = /(-!), «(1) = J=(1). 
Dadurch entsteht die Interpolationsformel 

m = j:(- 1) + *-±i [^(1) - H- 1)] + (*+;)(;-^) ^"(^,). 

Die Anwendung derselben gibt 

1 

//(Od< = ^(-l)+^(l)-|/"(T). 

— 1 

Daraus folgt die Näherungsformel 

b 



Jf{x)dx = ^^ \fia) + /•(&)] 



mit dem Fehler 

(&-a)\ r(£), 
1.23*12* 

Die rechte Seite der Näherungsformel ist der Flächeninhalt 
eines Trapezes. 

Nach Zerlegung des Integrals in n Summanden findet man die 
Naherungsformel 

b 

Jf{x)dx = ^ [(yo + yi) + (yi + y«) + •• •+ (y„-i + y„)]. 

a 

oder 



mit dem Fehler 



1 n* 1-2.2 12 



Wenn /^'(a?) dasselbe Vorzeichen hat für alle x zwischen a und 
hy so haben b^ und £ (in Nr. 22) verschiedene Vorzeichen. Deswegen 



ist das Integral j f(x)dx zwischen 
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. 2 2 2 2 

und 

enthalten. 

24. SüupsoiiBChe Formel. Wenn wir die Interpolations- 
formel (Nr. 15) 

S(t) = /(- 1) + *±i A H- 1) + ^^* A« j:(- 1) 

anwenden^ so wird 
1 

ß(f)dt = 2J:(- 1) + 2 AJ:(- 1) + i AV(- 1) - 4 J^^it). 

Es ist aber 
AJ:(- 1) = f(0) - J^(- 1), AV(- 1) = /(l) - 2/(0) + /(- 1), 

Nach einigen Reduktionen finden wir die Simpsonsclie Formel 



/' 



mdx^'-^[f{a) + ^f{^)+f(h)] 



mit dem Pehfer ^^__^, ^^^^ 



1-2-3-4-6 1-2-3-4 

Wenn f(x) eine ganze Punktion vom 3*®" oder niedrigeren Grade 
ist; so ist das Bestglied gleich Kuli und die Simpsonsche Formel 
ist ganz genau. Z. B.: 

(2a? + 3) (2a; - l)(2a? + 6)dx = 1 [- 3-5 + 5 -7] == y. 



Wenn aber f(x) eine andre Funktion ist, so kann man den 
Fehler verkleinem, nachdem man das gegebene Integral in n Sum- 
manden zerlegt hat. Man findet dann die Näherungsformel 



1 



a 

mit dem Fehler 



f{x)dx = -^ (i/o + 4y£ + 2^1 + 4^^ H f- 4y 2n-i + y«) 

2 2 2 



^___ 1 Q>-aY /"i^d) 



n* 1- 2-3. 4- 6 1-2-3.4 

In diesem Kapitel haben wir mit geringen Änderungen die 
Darstellung vom Herrn A. A. Markoff angegeben. (Differenzen- 
rechnung, S. 50 — 60.) 

Wir schließen damit die Anwendungen der Sätze über Differenzen 
und gehen zu den Summen über. 



Zweiter Teil. 

S n m m e n. 



Erstes Kapitel. 

TJnbestliiimte und bestimmte Summen. 

35. Definitionen, unter Summe von f{x) versteht man eine 
Funktion (p{x), deren Differenz gleich f{x) ist, also 

Diese Gleichung bedeutet, daß der Ausdruck qp(a; + A) — 9(^) 
80 transformiert werden kann, daß er in f{pc) übergeht. Dieselbe 
Relation wird auch bezeichnet durch 

,p{x) = 2f{x). 

Es ist vorausgesetzt, daß f{x) nicht in der Form ^{x + li)'—(p{x) 
gegeben ist Sonst wäre 

^{x) = 2[^{x + Ä) — qp(a?)] 

eine reine Identität, die keinen Nutzen haben könnte. 

Wir werden untersuchen, ob es möglich ist, daß eine andre 
Funktion ^(o;) dieselbe Eigenschaft 

H>{x + Ä) - H>{x) - fix) 

hat. Wenn wir ^(o?) — (p{x) mit F{x) bezeidmen, so ist 

i>{x) ^ ^{x) + F{x) 
und nach Nr. 3 

A^(a;) =» Ay(a;) + LI{x). 
Es ist also 

LF{x)^0 oder F{x + h)^F{x). 

Beide Funktionen ^(o;) und ^{x) unterscheiden sich also um 
eine periodische Funktion mit der Periode A. 
Ist 9{x) eine beliebige Funktion, fOr die 

0{x + Ä) = 9{x) oder A0(x) = 0, 
flo ist 

£^[(p(x) + «(a?)] = Aq>(x) + A 0(x) = f(x). 

SeliwAnoff, Dlfferenzeiuraohniuig. $ 
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Also haben alle Funktionen, deren Differenz gleich f{pc) ist^ die Form 

wo 9{x) eine willkürliche periodische Funktion mit der Periode h ist. 
Um die Aufgabe zu vereinfachen, werden wir von nun an voraus- 
setzen, daß X nur Werte der Reihe 

a, a + hy a + 2h, . . . 

annimmt. Dann ist 0(x) eine willkürliche Konstante, die wir mit C 
bezeichnen werden. Es ist also 

ufix) = (p(x) + a 

Dieser Ausdruck heißt unbestimmte Summe, weil der Wert 
von C ganz beliebig ist. 

26. Summatiou einflMher Funktionen. Aus jedem Satz 
über die Differenz (Nr. 2 — 6) können wir einen Satz über die Summe 
herleiten. 

Aus der Relation ^x = h folgt A^ » 1 und 



Weü 



2;i = f + a 



Ax(x — h)(x — 2h) ...(x — n — lÄ) = nhx(x — ä) . . . (a? — w — 2h) 
ist, so ist 

A :^x(x — h){x — 2h) . . . (a? — n — lÄ) = x(x — ä) . . . (a; — w — 2ä) 
und 
I!x(x — h){x — 2h) ..,(X'-n — 2h) = -j^x(x — h).,,(x — n — lh) + C. 

Damit man unter dem Summenzeichen n Faktoren hat, ist es 
zweckmäßig, n durch n + 1 zu ersetzen. Man findet 

2x(x — Ä) (rc — • 2h) ...(x — n — lh) = . j.±\k ^(?^ — ä) . . . (a: — nh) + C. 
Wir hatten die Formel (Nr. 2) 

= — nÄ • 



x(x + h)...ix + n-~lh) x(x + h).,.{x-\-nh) 

Im Nenner der rechten Seite sind mindestens zwei Faktoren^ 
weü w ^ 1 ist. 

Aus dieser Formel folgt 

x{x + h)...{x + nh) nh x(x+h) ,, .(x + n-lh) 

Die Funktion unter dem Summenzeichen hat im Nenner mindestens 
zwei Faktoren. 
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Wenn wir n — 1 statt n einsetzen^ bo entsteht die Formel 
« . 1 1 1 



x{x+h)..,(x + n-lh) (n-l)Ä x(x+h). ..(x+n-^h) 



+ a 



Hier ist w > 1. Dem Wert n = 1 entspricht eine besondere 
transzendente Fnnktion. die wir mit 27— bezeichnen. 

' X 

Die Relation 

Am* = Wf'^m^ — 1) oder A _ == m* 

gibt die Formel 

2Jw* = -7 — - + C. 

Es ist, • wie wir wissen (Nr. 6), 

A cos w = — 2 sm — sm ( w + -ö")* 
Wir setzen . ^ 

und wählen a mid /} so. daß 

wird. Ans der Bedingongsgleichnng 

ax + ß + 'Y='^ 
folgt ^ 

Deswegen ist 

A cos f a; — 2^] = — 2 sin 2^ • sin a?. 

Daraus folgt 

(-0 



und 



cos 



2sin — 
2 



= sina? 



cos I X 



(-1) 



2J sin a? = ^ — h (7. 

2sin^ 

Ganz analog findet man die Formel 

sin(a:--) 
2;cosrr = — ^—^ + C. 

2sin- 

27. Summation zusammengesetster Funktionen. Wenn 
A ein konstanter Faktor ist^ so ist' 

A-4.y(a:) =» -4.Ay(a;) 
^*^'' 2;[^Ag7(a;)] = Ä(p{x) + C. 

3* 
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Wenn wir A^(x) mit f(x) bezeichnen, so entsteht die Formel 

I]Äf(x)^AI]f(x). 

Die Konstante C ist schon im Ansdracke 2jf{x) enthalten. 
Wenn wir eine Reihe von Gleichungen haben 

A9?i(ir) « f^{x\ £^g>2(x) = f^{x), . . ., A^„,(x) = f^(x), 

so ist 

Daraas folgt 

^[/i(^) + f,(^)+"' + U(pd] - Viip) + f)P2(^) +• • • + f)Pm(:r) + G 
oder 

^[flip) + UiP) + • • • + /mW] = -2/i(a^) + ^AC^) + • • • + ^/m(4 

Mit Hilfe dieser Formel ist es leicht, ganze rationale Funktionen 
zu summieren. 
Es sei 



+ ^ 1.2.8...n 

Dann ist 

^/w-=^+A)nr"+^ — xT^TTh — +^ iTäTsÄ 

, , j (a; — a)(a? — a — Ä). ..(« — a — nÄ) 

■^ ^^ 1.2.8...(n + l).Ä ' 

Die Summe einer ganzen Funktion n*^ Grades ist also 
eine ganze Funktion vom (n + 1)*®^ Grade. 
Für Ä « 1 ist z. B. 

ic* = a? + a;(ic — 1), ar* = a? + 3a;(a; — 1) + ^{x — 1) (a; — 2) 

und folglich 

Zx^^G^\x{x--\)\\x{x^\){x-''2\ 

Zx^ == C J^^xix - 1) + a;(a; - l)(a; - 2) + ia;(a;- l)(a;-2)(ir-3). 
Es ist auch leicht, gebrochene Funktionen von der Form 

m 

(a: + »)(^ + ^ + ^) • • • (^ + * + **Ä) 
zu summieren, wo f{x) eine ganze Funktion ist. 

28. Partielle BmiunatioiL Aus der Formel 

t^\sp{x) • V'C^)] = SpC^) • A^(a;) + tp{x + ä) • Ag>(x) 

^^ U[g>(x) • A9(ir) + ^(a? + ä) • AyCa;)] = ^(a;) • ^(a:) + C 

oder 

2^g>(x) • A^(a?) + 2?9?(a; + ä) • /^^(x) = q){x) • ^(a;) + C. 
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Diese Gleichung liefert die Formel der partiellen Sammation 

Die willkürliche Eonstante erscheint, nachdem die Snmmation 
auf der rechten Seite ausgefQihrt ist. Mit Hilfe dieser Formel kann 
man das Produkt ^(x) • f(x) summieren, wenn man einen Wert von 
2f{x) kennt, und wenn die Operation 2J^(a; + Ä) • Ag?(ic) ausführbar 
ist, wo il}{x) = 2f{x) ist. 

Wir wollen z. B. 273*« rc für ä == 1 ausrechnen. 

Da ein Wert von ^3* gleich —3* ist, so ist 

2:3*. a; = 2:^; • A (| • 3*) = a; • |^ . 3*- 2;| . 3'+^ . Aa? 

Es wäre nicht^zweckmäßig, die Summation beim Faktor x aozufangen. 
Da Sx = -^-r — - ist, so wäre dann 

2:3*. a; = 2;3*- A ^^^""^^ = 3* • ?i?i:i) _ ^(Hii)^ ^3* 

m m a 

= |3*a;(a; - 1) - i:{x + l)a; • 3* 

und die gegebene Summe würde auf eine kompliziertere zurückgeführt. 
Noch ein Beispiel: 



2]x aiax = Ux-A 



-cos(a:--) 



2sm — 

2 J 



= X 



-co8(a:--) 



2 8m — 
2 



(2 8m-) 



2 sin — 2 sin — 

2 2 



29. Blgensoliaften der bestimmten BnmmeiL Wir haben 
schon gesehen (Nr. 25), daß, wenn 

g>(x + Ä) — q>(x) = /•(a;), ^(a? + ä) - f(x) = /"(a;) 

und a; eine Zahl der Reihe • 

a, a + Ä, a + 2Ä, ... 1 

ist, g>(x) und ^(a;) sich durch eine Eonstante unterscheiden 
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Subtrahiert man Yoneinander zwei Werte dieser Funktionen för 
x = a + nh und a? = a + wjä, so fallt C weg. Man findet die Gleichung 

g>(a + nh) — g>{a + mh) == ^(a + nh) — ip(a + mh). 

Der Ausdruck q)(a + nh) — y(a + mh) bleibt also ungeandert^ 
wenn man ^{x) durch einen andern Wert von 2]f(x) ersetzt. Dieser 
Ausdruck heißt bestimmte Summe und wird bezeichnet durch 

a -\-nh 

^,f(x) == g)(a 4- wÄ) — q)(a + mh). 

a-\-mh 

Die Zahlen a + mh und a + nh werden Grenzen dieser Summe 
genannt. 

Aus dieser Definition folgt* 

a +m h a+-(m+l)Ä 

^^Ca;) = 0, ^m = f{a + mh). 

Wenn c eine Zahl a + Jch ist, die zwischen a und h ^ a + nh 
liegt, so ist 

b 

^f(x) ^q>(b)- 9>(a) - [q,(b) - g>(e)] + [g,(e) - ,)(«)] 

a 

oder ft c ft 

Wenn man diesen Satz mehrmals anwendet, so findet man 

Daraus folgt eine wichtige Formel 
f{a) + f{a + Ä) + /•(« + 2ä; + ... + /^(6 - Ä) = <)p(6) - <)p(a). 

Hier ist f{pc) die gegebene Funktion, g?(a?) irgend eine Funktion, 
deren Differenz gleich f{x) ist, — ^ eine ganze positive Zahl. Diese 
Formel kann nur dann nützlich sein, wenn ^{x) nicht die Form hat 

f{a) + f{a + h) + f{a + 2h) + ^^' + fix^h). 

30. Berechnung einiger bestimmten Summen. 1) Es sei 

s == a 4- (a + Ä) + (a 4- 2ä) + h (» + w — lÄ). 

i. Wenn wir a + nh mit h bezeichnen, so wird 

b 

s^^x. 

a 
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Ein Wert der unbestimmten Summe Ux ist (Nr. 26) 

deswegen ist 

« = ÄtK6-Ä)-a(a-A)] 

oder nach einigen Reduktionen 

s = ^[a + (a + n- lA)]. 



18 



2) l-2-3 + 2-3-4+-+101112=2«('»-l)(«-2) 



= |-13121110. 

4 



88 



^ 1.8"^3-5"^5.7"^ "^^123 -^a:(aj + 2) 2 46 23 

n 
,4) l+2 + g«.+ ... + g«-l.=22'=^- 



5) 1 + cos Ä + cos 2Ä H h cos (w — 1)ä =^ cos a; 

o 

. 2n — 1, , . h 

sin — ^_Ä + sm- 

2sin — 
2 

Wenn wir von beiden Seiten — abziehen, so finden wir 

. 2n— 1, 



~ + COS Ä + cos 2Ä + • • • + COS (w — 1)Ä 



Bin 



2 



2 sin — 
2 



«Ä Ä 2n — 1, 

J*"^ cos — — cos Ä 

6) sin Ä + sin 2Ä H \- &in(n — l)h =^ sin a; = r^ 

2sin — 

2 

7) Wir wollen jetzt mit Anwendung von Formeln in Nr. 27 die 
Summen von Potenzen natürlicher Zahlen berechnen. 

X 

+ 1 + 2 + •■• + («- 1) =2* = ^^- 



X 



aj(aj-l)(2a?-l) 
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X 

08+ 18+ 2»+ ••• + (ic - 1)' --^^=^\x{x - 1) + x{x - 1) (a; - 2) 





— 1^12 



+ ia.(a.-l)(x-2)(a>-3) = p=^] 



Daraus folgt die bemerkenswerte Relation 
18+ 2»+ 3»+ ... + (a; - 1)»= [l + 2 + 3 + ..- + (a;- 1)]»* 



Wenn wir -t—tl luid t— tt-^ statt a? und a:* summieren, so finden wir 

1 • 2 1 • a • o 

X 

^'wIä ä* 1 «(« — 1) 



-i^iT2~" 1-2.8 2'l.2"^12^~" 1-2.3 ' 



^l-2*3 ■" 1.2-8. 4 2 *1.2.3'*'l2*1.2'~ 1.2-8.4' 



Im nächsten E^pitel werden wir die Funktion betrachten 



^1.2.8... (n-l) ~ ''"(^)' 

die, abgesehen von einem konstanten Faktor, von Jacob Bernoulli 
eingeführt ist. 

Zweites Kapitel. 
Die Jacob BemouUische Funktion. 

31. Bestimmimg der Koefflslenten. Nach dem in Nr. 27 
bewiesenen Satze ist die Jacob BernouUische Funktion 

X 

eine ganze Funktion n*^ Grades, die wir in der Form 

(2) y«(a')-Ar^^+^ i.2;'".(^_i) +A i.2.''.".(^_2) +--+A-ia? 

schreiben werden. Das konstante Glied fehlt, weil 9n(0) » ist 
Denn in Nr. 29 haben wir gesehen, daß die bestimmte Summe gleich 
Null ist, wenn ihre Grenzen einander gleich sind. Wir woUen die 
Koeffizienten ^, J^, ^, . . . , An—i bestimmen. 

♦ Vgl Encykl. der Math. Wias., Bd. I, p. 679. 
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Da 

x-{-l ^ X aj + l 4 

g>„(x + 1) ==^1.2... (n-l) ""^1.2... (n-1) +-2j 1.2...(n-l) 
ist; 80 besteht die Relation 

(3) g>„(x + 1) - 9„(a;) = ^.^'^'„_^y 

die nur f&r ganze positive x bewiesen ist. Die ganze Funktion 

Fix) = ip„(x + 1) - ip„(x) - 1.2!^, (^_i) 

vom nicht höheren als n*^ Grade verschwindet für jeden ganzen 
positiven Wert von x. Also ist F(x) identisch gleich Null und die 
Relation (3) besteht für jedes x. 

Wenn eine andre ganze Funktion ^(^) die Eigenschaften 

hat, so ist tl;(x) mit ^n(p^) identisch, weil beide Funktionen vom 
^ten Ghrade (Nr. 4) sind und außerdem die Gleichungen 

^(0) = ()p(0), ^(l) = ()p(l), i^{2)^^(2),.,.,tin)^q>(n) 

bestehen. Die Werte ^(l); 9? (2), . . . , g>{n) bestimmen sich sukzessive 
aus der Gleichung (3). 

Durch die Differentiation von (3) findet man 

y«('^+^)-y"(^)° i.2...(n-2) - 

Nach der Gleichung (2) ist g>n{0) == -4„_i. Deswegen ist 
Vn(p^) — ^—1 eiiie ganze Funktion, die für x==0 verschwindet und 

deren Differenz gleich 3-7-0 r ^9\ ^®*- Folglich ist für n> 1 

oder 

(4) <pn (x) = g?„_ 1 (ic) + ^1- 1. 

Aus dieser Relation schließen wir, daß 

ist; und daß 'die Koeffizienten J^, A^, . . . , An— 2 hier dieselben sind, 
wie in der Gleichung (2). 

Wenn wir in (3) a; = setzen, so finden wir g?»(l) = oder 
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Diese Relation gibt die Möglichkeit , Ät— i zu berechnen; sobald 
die Koeffizienten Äq, A^, . . ., An^i bekannt sind. 
Ans der Definition 

X 

folgt, daß \, 

% W 2 °° r2 ~ 2 * 

ist. Deswegen ist 

A=l, A 2' 

Die Gleichui^ (5) für n = 3, 4 und 5 gibt 

r-tT8 + i^ + ^ = 0' 



12. 3. 4-5 ' 12. 3-4 • 12-3 ' 12 

Daraus findet man 



12 ^ ' * 720 

In dieser Weise kann man beliebig yiele Koeffizienten berechnen. 

32. Bigensohaften der Koeffizienten. Wenn man in der 
Relation 

X durch — X ersetzt, so findet man 

(6) »»Si+lCl - «) - "Psi+l (- «) = i.2.8...2ft - 

Wenn man die Bezeichnung 

- Vst+iCl - «) = -^(«) 
einführt, so ist 

- »P^+xC- «) = - «Pit+itl - (^ + 1)] = Fio' + !)• 

Die Gleichung (6) nimmt die Form an 

^(^ + 1) - -^(^) = rrfrä*- 

Da aber jF(0) = ist, so ist nach dem in Nr. 31 bewiesenen 
Satze F(x) = 9'2jfc+i(^) ^^®^ 
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Vergleicht man die Ableitungen der beiden Seiten dieser Gleichung; 
so findet man mit Berücksichtigung der Gleichung (4) 

oder 

Für ifc = 1 ist die Ableitung von q>2k(p^) keine Jacob Bernoullische 
Punktion, weil die Definitionsgleichung (1) in Nr. 31 für w = 1 keinen 
Sinn hat. Aus der Differentiation von (8) folgt für 1c> 1 

oder nach (7) 

— -4.2t_i = -42*— 1. 

Es ist also 
(9) ^2ifc-i = für h>l. 

Wir haben schon in Nr. 31 gesehen, daß Ä^ nicht Null ist, und 
daß -4^ = isi Jetzt überzeugen wir uns, daß J^ = 0, -4^ = 0, . . . sind. 

33. Das Vorzelohen von 9>n(pc!) im Intervall (O, 1). 

Wenn man in (7) a; = ^ setzt, so findet man 



^»+1 (I) = 0- 



Wir wollen untersuchen, ob für ifc> 1 die Gleichung ^^g^+iC^)"^^ 
noch eine Wurzel im Intervall (0, 1) hat. 

Es sei vorausgesetzt, daß diese Gleichung außer und 1 noch 
die Wurzeln a^ und a^ besitzt, die den Ungleichheiten 

< ai < «2 < 1 

genügen. Hier darf a^ oder o^ gleich — genommen werden. Nach 
dem Rolleschen Satze hat 9?2jt i i(äj) = mindestens drei Wurzeln 
b^y &2 ^^^ ^3^ ^^ ^^^ Ungleichheiten genügen 

< 6i < »1 < ^2 ^ Ö&2 < 6^ < 1. 

Die Gleichung 9^2*4-i(^) '^ ^ ^** mindestens zwei Wurzeln c^ 
und C2, die zwischen und 1 liegen. 
Es ist aber 

weil ^2*— 1=0 ist. Die Gleichung 9^2t— 1(^) ^ ^ ^** außerdem die 
Wurzeln und 1. 

Aus der Voraussetzung, daß die Gleichung V^k+i^^) '^^ zwischen 
und 1 zwei Wurzeln hat, folgt also, daß die Gleichung 9^2* -iW'^^ 
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in demselben Interyall auch zwei Wurzeln liat. Wenn wir diesen 
Schluß mehrmals anwenden^ so finden wir, daß q>s(jx^) zwischen 
und 1 zwei Wurzeln hat. Das ist aber unmöglich, weil die Funktion 
3ten Qxades q)s(p^) i^cht für vier verschiedene Werte von x (0 und 1 
inbegriffen) verschwinden kann. 

Es ist also bewiesen, daß die Gleichung ?P2*+i(^)'^^ zwischen 

und 1 nur eine Wurzel x === — hat. 

Wir wollen jetzt die Gleichung (p^j^{oo) = untersuchen. Wenn 

diese Gleichung zwischen und 1 die Wurzel a hat, so folgt nach 
dem Bolleschen Satze, daß die Gleichung q>lj^ix) » mindestens 

zwei Wurzeln zwischen und 1 hat. Das ist aber unmöglich, weil 
9^2jfc(^) ^ ^'»jfc—iC^) ^®* ^^^ ^^ e\>en. bewiesen haben, daß ^^j^^^{p>) 
nicht ftir zwei Werte zwischen und 1 verschwinden kann. 

Es besteht also der Satz: Für alle Werte von x zwischen 
und 1 hat die Funktion (p^^(po) dasselbe Vorzeichen. 

Um dieses Vorzeichen zu bestimmen, betrachten wir die Gleichung 

9^2*+i(^) "" ^' ^® befriedigt wird f ür a; = und a? = g^- Nach dem 
Rolleschen Satze hat die Gleichung 

mindestens eine Wurzel x =^ a, die zwischen und ~ liegt. 
Aus der Relation 

folgt, daß Aii nicht gleich Null sein kann, weil ^^j^(x) nicht im 
Intervall (0, 1) verschwinden kann, und daß in demselben Intervall 
V2k(p^) ^^ entgegengesetzte Vorzeichen hat, als Äik* 
Es ist aber 

9^" W "" 12. ..2k "^ -^ 1.2...(2ä;-1) "^ ^ 1 • 2 . . . (2Ä-2) 
+ ^M.2...{2Ä-4)+***+^"-«r:^' 

Nach einem Satze [der Algebra hat g>2'i(x) ^ kleine positive 

Werte von x dasselbe Vorzeichen wie das Glied Ä^k^iT-Ti oder 

1 • 2 

wie Ä2k-i' 

Die Koeffizienten Äik—i ^uid Ä^k haben also verschiedene Vor- 
zeichen. Da aber A^ = ^> ist, so ist 

A^<0, A^>0, A^<0,... 

^"'^ ''* (- ly-'Ä^k > 0, 
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und far < a; < 1 besteht die Ungleicblieit 

(-l)*9'„(^)>0. 

Die Funktion Vjj • i(ic) wechselt ihr Vorzeichen im 

Interyall (0, 1), weil 00 = ^ eine einfache Wurzel der Gleichung 

g?j^,^(a?) = ist. Um diese Behauptung zu beweisen^ nehmen wir 
an^ daß 

sei. Die Funktion q)^j^(x) + -^a» verschwindet für a; == a und a? = -ö • 
Nach dem Rolleschen Satze soll 

verschwinden far einen Wert von x zwischen a und —f was un- 
möglich ist. 

34. BemonlUsohe Zahlen. Die Koeffizienten A^, Ä^, A^, . . . 
werden ofb durch die Zahlen B^^B^fB^,.,, ausgedrückt ^ die man 
Bernoullische Zahlen nennt. Es sind Zahlen^ die mit dem Vor- 
zeichen + oder — y als Koeffizient von x in der ersten Potenz, in 
den Ausdrücken für 

XXX 

Xi ^ } Xt ^y ^^ X y . , , 

auftreten. » 

Es ist 

X 

2^»* = 1 • 2t»(») = 1 • 2 A« - «* + y ' 


X 

2^=1-2 3.49?5(ar) = 1-2 3. 4^4ic+.-. 



X 

^a?«== 1.2- 3. 4- 5 -697(0;) = 1.2.3.4.5.6^^0; + ••• 


Es wird daher gesetzt 
1.2^ = jBi, 1.2.3.4ui4 = -J^, 12 3.4 5.6J[« = B3,... 

Allgemein ist 

(1^) ^" = i.8.8...2t - 

Die ersten Bernoullischen Zahlen sind 

i;,=._, i^=xj5» -^"42' ^i^W *5=66''*' 
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Die Tabelle der 62 ersten B findet man bei J. 0. Adams im 
Journal für Mathematik, Bd. 86 (1878), S. 269—272. Die ersten 
20 Zahlen sind im Lehrbuch von A.A. Markoff abgedruckt (S. 125). 

Auf die arithmetischen Eigenschaften der Bernoullischen Zahlen 
gehen wir hier nicht ein.* 

35. Bntwiokliing von ^^^aC^) ^ ®^® trigonometrisohe 
Beihe. Aus der Theorie der trigonometrischen Reihe (L. Dirichlet, 
Werke, Bd. I, S. 119 — 132) folgt, daß für 

0^a;<l 
die Entwicklung 

^^j^{x) == ö" ö^o + % ^^s ^^^ + ^ ^^^ ^^^ H h öf„ cos 2nicx H — 

+ \ sin 2%x + \ sin Axx -\ h 6„ sin 2nnx -\ 

besteht. Die Koeffizienten haben die Werte 

1 1 

an = 2 I q>2j^(x) cos 2njtX'dx, 6« = 2 / q>^j^(x) sin 2nxx • dx. 



Wenn wir x = l — y im Ausdruck für bn einsetzen, so wird 
wegen der Relation (8) in Nr. 32 



bn 2 / (p^j^(y) sin 2nity • dy. 





Folglich ist bn="-bn oder bn == 0. 
Der Wert für 



«0 = y SP2*(^)^' 



läßt sich aus.^er Relation 
bestimmen. Man findet 



> • > « 



.0 ... 

Um an zu berechnen, werden wir das Integral 

/* . /». . • . 

9^2* (^) ^^^ 2wÄir • da? durch / 9?2*-2(^) ^^^ 2w;ra5 • dx 

ausdrücken. Zu diesem Zweck gebrauchen wir die Formel 

* Diesem Gegenstand ist der dritte Abschnitt des Werkes von Herrn 
L. Saalschütz gewidmet (Vorlesungen über die Bernoullischen Zahlen . . . 
Berlin 1893). 
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(11) liptlt^dx =- q)tl>* — (p'tif + jil^ip'^dx, 

die eine Folge der zweimaligen partiellen Integration ist. 
Setzt man in (11) 

q> = q>^j^ {x), ^" = cos 2nnXy 

so findet man für h> 1 

f ' /\ , f sin2nÄaj 

Da aber 

1 



/< 



cos 2w5ra? • dx=^0 







ist, so wird 

1 1 



/ 9^2» (^) ^^^ 2wjti» • dic = ^^ — Tjj / 92t_2(a?) cos 2nxx • dx, 



Wenn man diese Formel {k — 1) mal anwendet, so erhält man 
1 1 

I 9^2* (^) ^^® 2w5ra? • drr = ■ 2k—2 I Vzi^) ^^^ 2n^x • dx. 



Es sei jetzt in (11) 

Daraus folgt 

9' = a;-2' ^'= ^^^' 

it ^ , C08 2n9rrr 



und 

1 



(2w«)« 



J ^^(x) COS 2n;riiJ • dx = ^^^r 





Man findet also 



(2wÄy 

Wenn man in der Entwicklunir 

• • • ' o 

/Pih(^) "? 2'-^o + ^1 ^08 27tx + öj COS 4;ric + • • • 
X gleich Null setzt, so findet man die bemerkenswerte Formel von Euler 

'^ ^ " (2«)" L 2" ^ 8" ^ 4" ^ J 
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Diese Formel wird später benutzt; am die Konvergenz einiger 
Reihen zu untersuchen. 

Aus den Werten fCLr A^ und A^ folgen die Formeln 

^ 2« ^ 3« ^ 4* ^ 6 

-*■ ^ 2* ^ 3* ^ 4* ^ ■ 90 

Dieses E^pitel habe ich unter dem Einfluß YonW.G.Imschenetzky 
(Easoner üniTersitäts- Nachrichten 1870 und französisch in J. Hoüel; 
Oours de calcul infinitesimal^ Paris 1888^ Bd. 1^ S. 476) und Herrn 
N. J. Sonin (Warschauer üniversitäts- Nachrichten 1888 und französisch 
im Journal für Mathematik, Bd. 116 (1896), S. 133) bearbeitet. 
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36. Bnlemohe Formel fOr gaase Panktloneii, Es sei 

Wenn man diese Funktion von bis x summiert, so findet man 






^f(x) ^PoVS^) + PlVn--M + • • • + Pn-jfPsC^) + 1>„-1^- 



Hier treten Bernoullische Funktionen ein. Wir wollen dieselben 
vollständig ausdrücken und die Glieder nach den Koeffizienten 
Ä^f Ä^y A^y . . . anordnen. 

In dieser Weise entsteht die Formel 

2f('')-Por7^+Pl l.2 .(n-l) + -' + ^n-,r2+Pn-l'^ 

^ ^ 

+ A(j>o i.2!'..(n-i) +i>i l.2!'.(n-2) + •^ + ^»-«^) 
+ ^ (i>o 1 . 2 !^. (n-2) +^x 1 . 2 '^. (n-8) + "• + Pn-B^) 



• • • • • • • 



+ An-iPaX 



Drittes Kapitel. EtQersclie Srnnmationsformel. 



49 



oder 



(1) 



X 



^ 



Das ist die Ealersche Sammationsformel. Sie ist genau für 
ganze Funktionen vom (n — 1)**^ oder niedrigeren Grade. 



37. Bestglied der Bnlemchen Formel. Wir wollen die 
Funktion 



(2) 



X j^ 

ifi^)-'^m- \MA<>^-A\f{^^-m'\-A^\f\x)-fm-■■■ 
. . • - A-.[/"(»-»>(a;)-/^(— »)(0)]- A-i[/'(»-')(*)-/'^"-'KO)] 



ausdrücken in dem Falle, daß f{x) keine ganze Funktion vom 
(n - 1)*~ Grade ist. 

Die Differ^iz von ^{^ ist 



(3) 



X 

A^(a;) = f{£) - A ff{x)dx - AAf(x) - A^dif'{x) 
J,_, A/'C— »)(a;) - A„-x6.f^«-*){x). 



Ans der Integralrechnnng ist bekannt, daß 



t^F{x) = Fix + 1) -F{x) = F\x) + ^F<\x) + ,^^F'"ix) + 



1 



U 



ist. 



Nachdem wir sukzessive 



X 

Fix) ^ff{x)dx, f^x\ f'(x), ..., f^-')(x) 



m = n 



n — 1, fi — 2, . . . , 



gesetzt haben^ so verwandelt sich (3) in 

Seliwanoff, Differenzenrechnimg. 



4 
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m -^^f'Qc) -^,nx) -... - ^-^/•(«-.(^) -/i^/^")(^+«)rf 







1 



Af(^)-^,n^)- ■ - i.2..Vl) ^^''""(^>- /^2':(n"-x/ ^"^(^+")^'* 





1 





1 
- A-i/"(— ') (x) - /a_i(1 - u)fi«)(x+u)du 



oder 



Af (a;) = - <p,(l)f'{x) - q>,(l)f"ix) - 

1 





Da die Werte 92(1); ^^^(l), ..., 9?„(1) gleich Null sind, so ist 

1 

Aflß(x) = — / VnO- — ^)f^''K^ + **)^«*- 


Nach der Formel (2) ist ^(0) = und folglich 

1 

(4) n,(x)=^-J^q>Sl-u)^f('')(x + u)du. 



Es wird hier nach x summiert und nach u integriert 
Wir werden n^2h setzen» Wegen der Relation 

ist dann 

1 

(5) i,{x) = - f9,,{^)du^r%x + u). 

Das ist das Bestglied toh Jacobi (Journal fOr Math., Bd. 12 (1834), 
p. 263 oder Werke, Bd. 6, p. 64). 
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Wenn man in (5) die Bernoullische Funktion ^^^.{yi) durch 
die trigonometrische Reihe (Nr. 35) ausdrückt^ so findet man das 
Bestglied von Poisson (Pariser Memoires, Bd. 6 (1826), p. 580). 
Wir werden die Jacobische Form gebrauchen. 

Die Funktion 9?j;^(w) behalt ihr Vorzeichen im Intervall zwischen 
imd 1. Deswegen folgt aus (5) 

1 



wo < < 1 isi 
Da aber 

80 ist 



9','i+i(«) = «PjiC«) + ^«»> 



/' 



1 



und folglich 



(7) 



(6) ti'(x) = Ä,t^f('^)(x + e). 



Wir haben also bewiesen, daß für jede Funktion f(x), die 
gewissen Stetigkeitsbedingungen genügt, die Formel besteht 

m + fii) + ^(2) + • . . + f(x - 1) 

X 

=jf{x)dx + A,[fix) - f(p)\ +Mf'i^) - rm +•■• 

+ ^,J, [/'(»-»)(«) - /'("-»)(0)] +^»_i [r("-»)(a;) - /'(»»-»)(0)] 
i+ ^2*[/'<"'(Ö) + /'<">(0 + 1) + • • + /^'"Kö + « - !)]• 
Für a; =- 1 ist 
1 

^0) ^Jf{x)dx + Mm - rm + A\fO) - rm + . • ■ 





(8) 



Die Glieder mit A^k^i fallen weg, sobald Ä; > 1 ist (Nr. 32, 
Form. 9). Für Ä; =- 1 ist 

4* 
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m+fii)+m+...+f(x-i) 

X 

(9) \ ^-Jmdx + Ä,[f(x) - f(Oy] 

{+\[r'(ß) + f'(ß + 1) + riß +2)+--+r(ß+x-i)] 

und 

1 

(10) ^(0) = ff{x)dx + Alfil) - f(0)] + A,f"(ß). 



38. Andre Form des Bestgliedes. Wenn die Ableitungen 
/^(")(0 und /•(«*+ >)(0 für Werte von t zwischen und a; ihre Vor- 
zeichen behalten und dieselben Vorzeichen habeU; kann das Best- 
glied (6) eine andre Form annehmen. 

In der Formel 

/ qxo^^du = (p(o' — (p*(o + I (Oip^'du 
setzen wir 

X 



Dann wird 

X 



X 

9" =^f^*'+'Kx + u), «, = 90„+,(«), 



und folglich ist 

-J^f'^K«' + «) • [9„(«) + A»]du 



» 



X 







. , Es besteht also die Relation 

1 







» 

1 1 

X 



° • 
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Da aber 
1 



/' 



f{n)(x + u)du - /'("-i)(a; + 1) - /•<"-i)(a?) = A/'("-i)(a;), 





so ist 



1 



fl»i')(x + u)du = f ("-«(«) - /"("-«(O) 

und folglich 

1 



-J*Vi,i^)ä»^f''"K^ + u) 





1 



Das erste und letzte Glied dieser Formel haben yerschiedene 
Vorzeichen^ weil^ nach dem in Nr. 33 bewiesenen Satze, T'siC^) ^^^ 
9jj , j(w) im Interyalle (0, 1) yerschiedene Vorzeichen haben und 
nach unsrer Voraussetziing f^^^\t) und f^^^^^\t) gleiche Vorzeichen 
haben. 

Die gefundene Relation hat die Form 

a = b + c, 
wo ac<0 ist. Aus der Gleichung 

a^= ab + ac 

folgt 

0<a^<ab, 0<|<1. 

Es ist also ^, 

wo < < 1 ist 

Nach dieser Bemerkung ist 

(11) fix) = . ÄittU^^^-Hx) - /'("-i)(0)]. 

Diese Formel kann statt (6) gebraucht werden, sobald f(t) gewissen 
Bedingungen genügt (Vgl. Markoff, Differenzenrechnung, S. 121.) 

39. Zweite Form der Enlerschen Formel. Wir werden 

yoraussetzen, daß für alle positiyen Werte yon t die Ableitungen 

r(i\ r\i), f^^i^), ' . • dieselben Vorzeichen behalten und daß aUe 

Ableitungen f(t), /"(t), f'"(t), . . . unendlich klein werden, sobald t 
unendlich groß wird. 
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In diesem Falle hat das Integral 

1 



A»(«)d«-2/'<"'(<+«) 



einen Sinn. Es bedeutet^ daß das Integral 

1 



ü 



bei unendlich wachsendem x sich einer bestimmten Gfrenze nähert. 

In der Tat, bei wachsendem x yergrößert sich der absolute 
Betrag des Integrals , weil neue Elemente mit demselben Vorzeichen 
zutreten. Sobald aber x hinreichend groß ist, so ist 



1 

ü 



<M„|{l + |/'("-«(0)|}, 



weil hier die Formel (11) anwendbar ist und nach der gemachten 

Voraussetzung \f^^^'^^^{o^\ beliebig klein sein kann. 

Nach bekanntem Prinzip hat das Integral eine bestimmte 

Gbrenze. 

Aus der Relation 

1 

« 

1 

= ^,»[/<"-^>(0) - /'<"-«(a;)] + /y,, . ,(«)d«2/<"+*(< + ~) 

folgt » 

1 

1 

f/ 



Da aber 
1 



q>,,(u)du2f(.*'')(t + u) und /9„+,(«)d«2/'(»+»>(< + u) 

*^ 



yerschiedene Vorzeichen haben, so ist nach der Bemerkung im 
Yorigen Paragraphen 
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00 



/i>„(«)dtt2'/'(»)(< + «) = 0,^»/-("-i)(O), 

wo < 01 < 1 ist. 

Um diese Formel zu Terallgemeinem, genügt es zu setzen 

wo V eine neue unabhängige Variable ist und x einen bestimmten 

Wert hat. 

Nach der bewiesenen Formel ist 
1 

r « 

oder 

1 

(12) ßMdu^f(^^)(t + u)^e, Ä,,f(^^-^^(x), 

* 

wo < Öa < 1 ist. 

Wenn man in der Eni ersehen Formel das Bestglied 

1 

^ 

durch 

1 1 



ersetzt; so findet man 

^m = lfii)ät +Ct + Ä^f{x) + Ä^f{x) + ■■■ 

*/ 



» 



* 

wo 

a = -^r(0)-^r(0) — •• 

1 

/ 

Jetzt wollen wir beweisen^ daß die Eonstante Ck von h un- 
abMngig ist. 
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Für Ä = 1 ist 

Z 



Mit Benutzung Ton (12) findet man 

X 
X j% 

und , ^ 

Der Yei^leich der beiden Formeln gibt die B«Iation 
Ci + A,r(x) + -- + ^,t_xf(»-»)(a;) + ^,t /'(«*-« (a;) 

= C, + 0,^A«), 

die für jeden ganzzahligen positiven Wert von x besteht. 

Bei unendlich wachsendem x sind alle Ableitungen Ton f(xy 
imendlich klein und folglich ist Gk » C^. 

Es ist also die Formel bewiesen 

X 

(13) ^m-G^jmdt^AJ{x)^A,f{x)^- 

wo < < 1 ist. 

40. Der allgemeinere Fall der Bnlersohen FonueL Bis 

jetzt haben wir die bestimmte Summe Ton der Form 

m+m+f(2)+...+f(x-i) 

betrachtet. Nun wollen wir zur Summe 

F{a) + F(a + ä) + l^(a + 2ä) + • • • + F{b - h), 

wo b ^ a + xh ist, übergehen. 

^ '*' F(a + th) = ^(0; 

dann ist 

F(p - A) = fix - 1) 
und folglich 

F(a) + F(a + h) +'■■ + F(b-K) -fiO) +fO.) +-- + fix- 1), 
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Wenn wir berücksichtigen, daß 

f{fn)(f) = A"»JF'<'»)(a + th), 
/•(»)(«) = Ä«».F('»)(6), /•(»)(0) = Ä'»JF'(«)(a), 

80 nehmen die Formeln (7), (8), (9) nnd (10) die Formen an 
•F(a) + F(a + h) + F(a + 2h) + --' + F(h — K) 
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(14) 



^rFiu)du+Äi[F{b)-F(ay]+A,h[F'{h)-F'(a)] + 



(15) 



+ Än-i A»*-» [J'("-»)(&) - J'<»*-»)(a)] 

+ ^„Ä"[l?'(»»)(a+eA)+J'<")(a+l+0A)+- • ^-.^»»(fe-A+eÄ)], 

a+k 

F{a)=^rFiu)du+AilF{a+h)-F{a)]+A^hlF'(a+h)-F'{ay]+- ■ ■ 

+^,i_iÄ»*-*[J'(«t-«)(a+A)-J'("-»)(o)]+^,iA"l?'{«*)(a+0Ä), 
'F(a) + F(a + h) + F(a + 2h) + "- + F(p-h) 

b 



(16) 



= ifF(u)du +A, [F(h) - F(a)] 



nnd 



+ J,A» [^"(a+eÄ) +^"(a+ 1+ 0A) + • • • +.F"(6-Ä+eÄ)] 
(17) ^(a) - ^ JF(u)du + ^1 [.F(a + A) - J'(o)] + J,Ä» J"'(a + 0Ä). 

a 

In allen diesen Formeln ist h ^ a + xh und rr eine ganze 
positive Zahl. 

Wenn die Ableitungen JP"(0, F^(f)y F^^)(t), . . . fOr < ^ a die- 
selben Vorzeichen behalten und die Ableitmigen F\t), F*\t), F'"(t)y.. . 
unendlich klein werden^ sobald t unendlich groß wird, dajin besteht 
die Formel 



(18) 



{F(a) + F(a + h) + F{^a + 2h) + --- + F(b-h) 

b 

- (7 + ^ CF(u)du +^1 FQ,) + J, Ä 1?"(5) + • . • 
+ Au-i Ä«*-« J'(«*-«)(6) + 0^,t A«*-i2^("-i)(6), 



wo die Konstante C Ton h und % unabhängig ist. (Vgl. Markoff, 
Differenzenrechnnng, S. 131 — 133.) 
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Viertes EapiteL 
Anwendungen der Enlerschen FormeL 

41. Sntwicklnng von —^ — in eine Beihe. Setzt man in 

e — 1 

der Formel (8) (Nr. 37) 

80 wird 

1 



/ 



e"dx = -j^, /•(»)(!) = f-e«, ^(«)(0) = <«, /-(»«(x) = <«e' 



und folglich 

l = ^+^i(e'-l)+^<(e'-l) + -- 
+ A,i-it"-'(e< - 1) +^„_i<"-»(e'- 1) + ^,t<«*e»'. 

Nach der Formel (12) (Nr. 35) ist das Resi^lied 

Bei wachsendem Je bleibt 

endlich, weil 

oder t 

l<Jf*< i^- 



1- 



gSA-l 



ist. Daraus folgt, daß lim Mi == 1 ist. 

Es ist also für |^|<2^ das Restglied dem absoluten Betrage 
nach beliebig kleiu fiir hinreichend große Werte Yon Je 

Folglich besteht für | ^ | < 2« die konvergente Reihe 



e*-l 



• • m 



1 = ^ + A(e'- 1) +AKe' - 1) +^«<»(c' - 1) + 
oder 
(1) -^^=j+A^ + A,t + Ä^i»+-- 

e — 1 * 

42. Entwicklung von ctg x und tg a? In eine Beihe. 

In der Formel (8) (Nr. 37) setzen wir 

Dann wird ^^^^ = '''' ('^> 



1 



/ 



(iOB{tx)dx = ^^, f'(t)='-t6m(tx), f"'(t) = + i^Bm(tx), 
. . ., fi.i>ö(x) = (- l)*-i«» cos (tx), 
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und folglich ist 

1 =« -^ + A^ (cos ^ — 1) — A^t sin ^ + Aji^ sin ^ 

+ .l2i-2(- l)*-i^2*-8 sin t + ^2,(- l)*-i^«* cos (00. 
Da aber 

•< A / X i\ 1 + COS t ^ t 

1 — uii(C0S * — 1) = —^ = COS^ -; 

SO findet man nach der Division durch 

sin ^ = 2 sin — cos — 
die Formel 

Ebenso^ wie in Nr. 41^ kann man beweisen^ daß das Bestglied 
für \t\ <.2jc und für hinreichend große Je beliebig klein wird. 

Wenn man x statt — einsetzt^ so entsteht die Entwicklimg 

(2) ctg ii; == - - 2^A^x + 2^A^x^ - 2^A^x^ + • • •; 

die für | o? | < st konvergiert. 
Mit Hilfe der Relation 

tg 0? = ctg X ■— 2 ctg 2x 
findet man die Reihe 

(3) tg ic = 22 (2>- 1)^0? - 2\2^-^l)A^x^ + 2^(2^-l)A^x^ ; 

wenn | a; 1 < -r^ ist. 

43. Stlrlingsohe Beihe. Wir wenden die Formel (18) (Nr. 40) 
an und setzen 

F(x)^logx, a = \y Ä = l, 'b = x. 
Man findet dann 

/ log xdx = X log X — X, F\x) = a?""^, F^^(x) = — or-^, • • • 

..., i^M(a;) = (-l)"»-U.2.3...(m-l)^'» 
imd folglich 

log [1 • 2 • 3 . . . (a? — 1)] = C + i» log X-- X +Ailogx 
+ J^izr-i + ^4l.2izr-8 + -.- + ^2A_.,1.2.3...(2Ä-4)ii;-<«*-8) 

+ 0.l2*l-2-3-..(2Ä-2)ar-(2*~i). 
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Indem man zu beiden Seiten log x addiert und Ä^ durch 
seinen Wert ( — — j ersetzt, so erhalt man 

log (l '2 '3 . . . x) == C + (x + ^"jlog X -' X + A^x-^ 

(4) |+^J.2a?-» + .-- + ^,ife_,1.2.3...(2*-4)a?-("-») 

+ 0^2* 1 • 2 . 3 ... (2* - 2)ip-("-i). 

Hier ist Z; > 1 Torausgesetzt Für Jc==l besteht die Formel 

(5) log (1^2 'Z ... x) =- C + (x + ^"^log X - X + e^A^x-K 

Es sind 6 und di Größen zwischen und 1. 
Es bleibt noch die Eonstante C zu bestimmen. Zu diesem 
Zweck wenden wir die Wallissche Formel an: 

— » lim tf (n) für lim n = oo, 
wo 

^/x 2*2 4-4 6-6 2n*2n 



1.8 8.6 6-7 (2n-l)(2n+l) 

ist. Mittels einfEtcher Umformungen erhalten wir 

^^^^ Ll-8-6...(2n-l)J 2n + l 

_ r2'2.4.4.6.6...2n.2n "[» 1 
""L l-2-8.4...(2n-l)2n J '2n + l 

^ (2.4-6... 2n)^ 1 _ 2*** . (1 • 2 . 8 . . . n)^ 1 

^ (1.2.8.4...2n)* ' 2n + l^ (1 • 2 . 8. . .2n)«.2n \ , 1 " 

^ + 2^ 

Daraus folgt 
log *(n) = 4n log 2 + 4 log (1 • 2 . 3 . . . n) 

- 2 log (1 .2.3 ... 2n) - log (2n) - log (l + ^) 

oder mit Benutzimg Ton (5) 

log *(n) = 4n log 2 + 4(7 + (4n + 2) log n - 4n + 461 J, ~ 

- log 2 - logn -log(l + ^) 

- (4n + l)log 2 - 2C- (4n + l)'logn + 4n - 6, J,^ 

• * » 

21og2 + 2(7+4e,Ai-log(l + ^)-e,^~ 
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Wenn wir jetzt n unendlich wachsen lassen, so wird 
lim log *(n) = log |^ lim [461 ^ -^J = 0, 

1^1^8(l + Ä)=0, lim[0,^i] = O 

und folglich 

log|-~21ög2 + 2a 

Daraus findet man 

= log>/2är. 

Es bestehen also die Formeln 

log ( 1 . 2 . 3 . . . rc) = log yW^c + (x + ^\log X — X + Ä^x-^ 

(6) \+Ä^l'2x'-^ + " +^2t_2l.2.3...(2*-4>-("-») 

+ 6^2*1 • 2 . 3 . . . (2Ä; - 2)a;-(«*-i) 
und 

(7) log (1 . 2 . 3 . . . rc) = log >/2i + (x + ^^ogx - x+ d^A^ x^K 

Da < 01 < 1 und J, = ^ ist, so folgt aus (7) 

log y2n; + (rc + i) log ir — rc < log (1 . 2 . 3 . . . rc) < log >/2^ 

oder • • ... . . ^ 

(8) y2ixx • ÄJ* • C-* < 1 • 2 • 3 . . . a? < y2n;x • a?* • c-* • e^. 

Für große Werte Yon x kann die Naherungsformel 

(9) 1 • 2 • 3 . . . rc = y27cx - af - e-' 
gebraucht werden*. 

44. Vntersaohnng des Bes^liedes der StlrlingBohen 
Beihe. Die Gleichimg (6) ist die Stirlingsche Reihe, deren Eon- 
yergenz wir jetzt untersuchen * wollen. Zu diesem Zwecke muß man 
sehen, ob das Restglied bei unendlich wachsendem Je unendlich 
klein wird. 

Nach der Formel (12) (Nr. 35) ist 

0^2* -1 • 2 . 3 ... {2h - 2)ic-("~i) 

"^ \ ^^ ^"^ \^%i] 2«aj ««a? L ^ 2" 8** 4'* J 
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In Nr. 41 ist bewiesen; daß bei wachsendem k sich die Reihe 

^ ' 2^* 8** 4** 

der Grenze Eins nähert. Für jeden Wert von x wird aber das Produkt 

12 8 2Ä;-2 



• • • 



%%x 2%x ^%x ^%x 

beliebig groß für hinreichend große Werte von Tc, 

Daher ist die Stirlingsche Reihe divergent. Für große Werte 
von X und kleine Werte von Tc ist das Restglied dem absoluten 
Betrage nach sehr klein. In diesem Falle gibt die Stirlingsche 
Reihe eine gute Annäherung für 

log (1 • 2 . 3 . . . :r). 

Z. B. für a; = 100 und i = 3 ist der Fehler < — ^- (Markoff, 
Differenzenrechnung, S. 137.) 

Für einen bestimmten Wert von x ist der Fehler dem absoluten 
Betrage nach kleiner als 

1 .1,* 1 . 1 . 2 . 3 ... (2* - 2)a;-(»*-i). 

Die beste Annäherung bekommt man für einen solchen Wert 
von ICf für den der Ausdruck 

Uk = 2X' [j^j • 2^ • 2^ 2«S" L "^ ?* "^ i" "^ * * J 

den kleinsten Wert annimmt, weil der absolute Betrag des Restgliedes 
kleiner als U]f ist. 
Für Ä ^ 6 ist 

' + — + — + — + •••< ' 



Deswegen kann man angenähert setzen 

1 \8 2 8 2Ä;-2 



• • • 



Uh = 2x • ( 1 • • — 

X^TCxJ 2^x 2^x 27CX 

Da aber angenähert 

— 9 / 1 y 2 ^ 2fe~2 ^ 2Ä;~1 ^ 2Ä; 
SO ist 

Wt 2«a5 2nx \xx/ 
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Die absoluten Beträge der Glieder der Stirlingschen Reihe 
nehmen also ab; solange Je <.7CX bleibt. Die beste Annäherung erfolgt 
also, wenn Je nahe zu jra? + 1 liegt.* 



45. Angenäherte Berechnung der Summe 

fftr groClEie Werte von x. In der Formel (18) (Nr. 40) setzen wir 

F(x) = —} a = l, Ji = l, h=^x. 

Dann finden wir 

j^ = \ogx, -F('«)(a?) = (-l)'»1.2.3...w^-(«+i) 
und folglich 



(10) 



1 + ^ + 1 + - + ^ = C + log 0^-1 ^-^o;-^ 

-.lj.2.3a;-* 42*-2l-2.3...(2Ä-3)a:~("-») 

-e.l2ifcl-2.3...(2Ä-l)a?-". 



Die Konstante C läßt sich nicht durch bekannte Eonstanten aus- 
drücken^ sie läßt sich nur angenähert berechnen. 

Wenn man x =» 10^ Jc^l setzt imd das Bestglied yemachlässigt; 
so findet man 

0« 0,57721 56649 01532... 

Das ist die Eulersche Konstante. 

Ebenso wie in Nr. 43 kann man zeigen^ daß die Reihe (10) für 
jeden Wert von x divergiert. Deswegen darf zur Berechnung der 
betrachteten Summe nur eine geringe Anzahl der Glieder der Reihe 
genommen werden. 

In die Untersuchung^ wie groß diese Anzahl sein kann, werden 
wir uns nicht einlassen. 



♦ Vgl. EncyU. der Math. Wias., Bd.I, p. 108. 
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Erstes Kapitel. 

Allgemeines über Differenzengleiclmiigeii. 

46. Zwei Formen der Differensenglelchimgen. Eine 
Gleichnng von der Form 

(1) *(^; y^j ^y^y ^*y«; •••, ^""yj = o 

heißt Differenzengleichung. Hier ist ^ eine gegebene Funktion, 
y eine unbekannte Funktion. 

Bei Anwendung der Formeln 

^% = y.+2 - ^y.+i + y.^ 



^"y. - y.+fn - ^y^+m-i + ^^^rV^ y^+m-» — + (- ^Tv^ 

nimmt die Gleichung (1) die andre Form an 

(2) tl;(x, y^, y^_^^, y^+2, ..., y^+J = 0. 

Ist die Gleichung in der Form (2) gegeben^ so kann man sie in 
die Form (1) überfähren. Es genügt, die Formehi 

y.+i-y.+ ^y,y 



y.+n» = y. + ^^y. + ^^^tV^ ^'y« ^ 

anzuwenden. 

Von diesen beiden Formen (1) und (2) ist (2) die bequemere. 
Es ist einfacher, eine gegebene Funktion statt y^ in die Gleichung (2), 
als in (1) einzusetzen. 
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Wenn bei der Transformation (1) in (2) z. B. die Größen y^ 
und y^ , ^ wegfallen^ so wird die Gleiclmng (2) dadurch vereinfacht, 
dafi man y^,«=;ef setzt. 

Es sei die Gleichung 

^'y, + ^% - Ay, - y, = 
gegeben. Mit Anwendung der Formeln 

^\= y.+»- 3y»+»+ 3y,+i- y,, 

geht diese Gleichung in 

über. Setzt man yg.t^'^^xf ^^ entsteht die einfachere Gleichung 

Aus diesen Gründen werden wir YoraussetzeU; daß die gegebene 
Gleichimg die Form (2) hat und daß in ihr y^ enthalten ist. Wenn 
außer y^ noch y^ , ^, aber kein y mit größerem Index vorkommt, so 
werden wir (2) eine Gleichung w**' Ordnung nennen. Die eben 
betrachtete Gleichung 

^%+^%-^yx-yx^^ 

ist also eine Gleichimg 1*®' Ordnung. 

• • 

47. Partikulftre und allgemeine Lösungen. Wir werden 
voraussetzen, daß die vorgelegte Differenzengleichung n^^' Ordnung 
die Form hat 

Die Funktion f soll für alle betrachteten Werte von a?, y^, y^j+i? 
Vx+v'"^ Vx+n—i ®^ö^ endlichen und bestimmten Wert erhalten. 

Jeder Ausdruck y^=^(x), der dieser Gleichimg genügt, heißt 
partikuläre Lösung von (3). 

Es gibt imendlich viele partiku^e Lösungen, um dies zu zeigen, 
wollen wir den Größen y^, y^, y^,", V ^i willkürliche Werte 
beilegen. 

Aus (3) ergibt sich für a? = 0, 1, 2, 3, ... das System der 
Gleichungen 

Seliwanoff, Differenzenrechnung. 5 
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Vn-^fi^f Voy Vu y%y''y »„-1)7 

Vn+i^fQ^y vty y^y y^y-y y^y . 

»«+,"- /"(S, %; »4, »67..., y„+,)' 

aus welchem y^, y^^^, »«+2^ »n+s^ • • • s^J^össi^ö berechnet werden 

können. Es hat also y^ einen ganz bestimmten Wert fßr ganze 

positive Xf sobald y^^ y^j y^y.-j y^_^ gegeben sind. Man kann also 

schreiben 

(5) y.= F{Xy y^y y^, y„ . . ., y^__^). 

Die Wahl der Anfangswerte y^, y^, %,..., y„_i ist nur einer 
Bedingong unterworfen: es sollen alle Gleichungen des Systems (4) 
einen Sinn haben. 

Es ist möglich, daB die Funktion 

f(!^> y»^ y.+v yx+v'-y y.+«-i) 

einen bestimmten Wert erhält für a;^a und fär beliebige Werte 
▼on y^y y^+i, y,+„..., y^+^-r 
Dann besteht das System 

ya+n^f{<^y ya> Va+V »a+J^"-^ »a+n-l)^ 

»«+, + !--/'(«+ 1; »a+17 ya + V Va^V ' ' f ya + n)> 
ya^n + 2^f{^ + ^> »a+2> ^a+B' »«+4^ ' ' • > »a+n+O^ 



(6) 



Daraus folgt die Relation 

&Lr x'^a. Hier ist die Differenz x — a gleich einer ganzen posi- 
tiven Zahl. 

In den Formeln (5) und (7) sind y^, y^, y^, ..., y^^i oder 

ya> Va+xy Va+v '> »aH-n-i wiUkürUche Konstanten. 
Jede Lösung der Gleichung (3) von der Form 

(8) y,^t{Xy C^y C„ C„ ..., OJ, 

die n willkürliche Eonstanten enthalt; heißt allgemeine Lösung, 
wenn es möglich ist, die Gleichung (8) in die Form (5) oder (7) 
zu bringen. 
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Zu diesem Zweck muß man aus dem System 



yi = ^(l, Ci, Ca, ..., (7„), 



(9) oder (10) 






Ci, Cg, . . ., c^, 



^y«H-n_i=*(ö^+W-l,Ci, c'a, ..., cj 



die Konstanten Ci, Cg, . . • , C^ durch y^j, Pv V^f -- > Vn-i ^^^^ durch 
»a^ ^a+i' ya+2> ••'^ ^«+»-1 ausdrücken. 

Man nennt die willkürlichen Eonstanten C^i, Cg, . . ., Cn wesent- 
lich, wenn das System (9) oder (10) in bezug auf diese Eonstanten 
lösbar ist. 

Allgemeine Lösung von (3) ist also eine Lösung, die n 
wesentliche Eonstanten enthält. 

Aus einer allgemeinen Lösung lassen sich viele andere ableiten. 
Es genügt in (8) 

(Ä== 1, 2, . . ., n) 
zu setzen. Die Eoeffizienten a^^ sind nur einer Bedingung unter- 



worfen, daß die Determinante 



a 



kh 



(fc = 1, 2, . . . , n; Ä = 1, 2, ...,») 



nicht verschwindet. Li diesem Falle lassen sich c^, o^, . . . , o„ durch 
(^v Cgi • • • ; Cn ausdrücken. 

48. Eigensohafben der linearen Dlfferenzengleiohnngen. 

Eine Gleichung von der Form 

(11) y.+„+i>iy,+„-i+Ay.+n-2+-+i>n-iy.+i+i^«y.= ö, 

heißt lineare Gleichung. Die Eoeffizienten j)^, p^, . . ., i>„_i> i>„ 
und das letzte Glied Q^ sind gegebene Funktionen von x. Wenn 
wir Q^ durch Null ersetzen, so entsteht die entsprechende 
homogene Gleichung 

Die Gleichung (11) heißt eine vollständige, wenn Q^ nicht 

gleich Null ist. 

Es sei Tx eine partikuläre Lösung von (11). Wenn man in 
die linke Seite dieser Gleichung 

(13) y.= ^.+^x 

einsetzt^ so findet man 

5* 
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Die erste Zeile ist nach der Voraussetzung gleich Q^. Daraus 
folgt; daß der Ausdruck (13) der vollständigen Gleichung genügt, 
wenn 0^ euie Lösung der entsprechenden homogenen Gleichung (12) ist. 

allgemeine Lösung von (12) ist, dann ist 

V:,^ ^«,+ Vi^7 ^p Ci» . • •; C'*) 

allgemeine Lösung der vollständigen Gleichung (11). 

Lineare homogene Gleichungen besitzen folgende bemerkenswerte 
Eigenschaften: 

1. Wenn z^ eine Lösung von (12) ist, so ist auch Cz^ 
eine Lösung. Dies folgt aus der Relation 

= C'C^.+n + i>l^«+«-l + A Vn-2 + • • + Pn^.y 

Hier ist C eine willkürliche Eonstante. 

2. Wenn z^'^ und zj^^ Lösungen von (12) sind, so ist auch 
zjp + z^^ eine Lösung, weil 

+ pM' + ^'^O = {'fU + Pi'^U-^ + Ä^Ä„_, + • • • + i'-^n 

ist. 

Aus diesen Sätzen folgt der wichtige Satz: 

3. Sind z^^^j z^^'^. .... z^^^ verschiedene Lösungen der 
linearen homogenen Gleichung, so ist auch 

(14) y,= C,0^) + 0,;.W + • • • + Cy«) 

eine Lösung. Hier sind C^ C^, . ^ ., C^ willkürliche Eonstanten. 
Wenn dieselben wesentlich verschieden sind, dann ist der Ausdruck 
(14) allgemeine Lösung von (12). 

Li Nr. 47 haben wir gesehen, daß eine Differenzengleichung eine 
Funktion definiert, wenn die Anfangswerte gegeben sind. XTm n 
partikuläre Losungen zu bilden, genügt es, n Systeme von Anfangs- 
werten anzugeben 
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(15) 



^'O f 1 ^ ^2 ' ' n — V 

ein) g{n) An) Mn) 



Wenn die Determinante^ die aus diesen Elementen gebildet ist, 
nicht yerschwindet; dann sind im Ausdruck (14) die Konstanten 
wesenÜich yerscliieden und die Formel (14) ist wirklich allgemeine 
Lösung der Gleichung (12). 

Wir bekommen also den Satz: 

Die allgemeine Lösung einer linearen homogenen 
Gleichung ist eine lineare homogene Funktion der will- 
kürlichen Eonstanten. 

Es war bei dem Beweise dieses Satzes vorausgesetzt, daß die 
Koeffizienten der gegebenen homogenen Gleichung einen Sinn haben 
für ganze positive x. 

Wenn aber wir nur behaupten können, daß i)i, i>2; • • •? i>„ ^^'^' 
liehe und bestimmte Werte haben für rr ^ a, dann ist das System (15) 
durch das folgende 



(16) 



a ' a+1' f a-fn — V 

a ' a+1' ' a-f-n — V 









' } 



^in) 



ZU 



ZU ersetzen. 

Nach diesen allgemeinen Betrachtungen gehen wir jetzt 
speziellen Gleichungen über. 

Zweites Kapitel. 
Lineare Differenzengleichnngen erster Ordnnng. 

49. Homogene Olelohnng. ITm die allgemeine Lösung der 
Gleichung 

zu bilden, geben wir der Größe yQ einen beliebigen Wert. Aus 
(1) folgt dann 



y«"~"-^-^l-^ • • • -^a;— 1^0- 
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Weil j/q willkürlich ist; kann man schreiben 
(2) y^^C'A^'Ä^'A^ ... Ä^_^, 

wo C die willkürliche Eonstante ist. Das ist die gesuchte Lösung. 

Wir werden voraussetzen, daß in (2) keine der Gh-ößen A^, Ä^, 
A^, ..., Äx^i gleich Null isi Sonst hätten wir die evidente 
Lösung y^ = 0, die kein Interesse bietet. 

Wenn für alle x'^a die Funktion Äx einen endlichen, be- 
stimmten und von Null verschiedenen Wert hat, dann ist die 
allgemeine Lösung unserer Gleichung 

(3) y^=-G'Äa'Äa+l'AaJ^2 ... 4r-l- 

50. Vollst&ndlge Olelohong. um die Gleichung 

(4) y.^i-Ay.-^. 

aufzulösen, setzen wir 

(5) y^^-UxVx, 

wo Ux und Vx neue unbekannte Funktionen sind; eine derselben 
wählen wir so, daß die Gleichung für die andere Funktion möglich 
einfach sei. 
Da 

Vx+l^Vx + £^Vx 

ist, so nimmt die Gleichung 

Ux^l • Vx^l — AxÜxVx = Bx 

die Form an 
oder 

(6) (Ux^i — ÄxUx)Vx + Ux+i • £^Vx = Bx. 

Wenn Ux eine partikuläre Lösung der entsprechenden homogenen 
Gleichung (1) ist, dann erhalten wir die Gleichung 

Ux^l'^Vx = Bxf 

die Vx bestimmt. Es ist 
und folglich 

Aus der Betrachtung in Nr. 49 folgt, daß man 

Ux = A^Ä^A^ . . . Ax—i oder Ux = AaAa^i . . . Ax—i 

setzen kann. Jn diesen partikularen Lösungen ist ti^ =» 1 oder Ua =° 1. 
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51. Beispiele. 1. Es sei die Gleichung 
gegeben. Hier ist -4« = 3, Bx^ »*, 



Die Formel 



^-=2'5ä:i + ^ 



hat Terschiedene Werte, je nachdem a = 3 oder a ^ 3. 
Wenn a von 3 verschieden ist, dann ist 

8 
und folglich 



a' 



Wenn aber a =-= 3 ist, so ist 

und 1 

y,==C.3* + i^.3-. 

2. Wenn wir die Gleichung 



nehmen, so ist 



, 3a?+l X 






3aj + 7 ^' (3a; + 4)(8aj+7)' 

7 10 3« — 5 3aj — 2 



7 10 13 16 3aj+l 3« + 4 



= (_l)x. !•* 






(3a!+l)(8a!+4) 
i(-l)«.a;, 



2;(- 1)' . a; = 2;ar . A (^ = - \x{- 1)»- -^2;(- 1)' 
Folglich ist die allgemeine Losung dieser Gleichung 

^«""("* ^)*'(8aj+l)(3aj + 4)' "^ U^""!/ ' (3a;+l)(8aj + 4) 
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3. Als letztes Beispiel wollen wir noch die Gleichung 

^ 2x + l Zx + 1 1 

auflösen. Es ist 

7 9 11 2x + 6 
^ 13 5 2a;- 1 

14 7 3a:-2 (2a; + 1) (2a; +3) (2« + 6) 



4 7 10 3a; + l 15(3a; + l) 

3__ 16 

^a.4.1 (2a; + 3)(2a; + 5)(2a;+7)' 

^(2a; + 3)(2a; + 5)(2a;+7) ^ 8"^' 
1 1 



•; 



(•+i)('+i) 



K)KI)('+I) 



10/ I ^\/ 1 ö\ ' 4 (2a; + 3) (2a; + 6) 



^ (2a;+l)(2a;+3)(2a; + 6) ^_1 2a; + l 
^x"" 3a; + l 4 * 3a;-fl' 

(Vgl. Markoff, A.A. Diflferenzenrechnung, S. 156 — 158.) 
Jetzt wollen wir die Anwendungen der Gleichungen erster 
Ordnung zeigen. Die Koef&ienten verschiedener Entwicklungen 
genügen Differenzengleichungen. Die Auflösung derselben führt 
zur Bestimmung der Koeffizienten. ' Als Beispiel nehmen wir die 
binomische Entwicklung und die Entwicklung von cos xt nach Potenzen 
von cos t 

52. Die binomisohe Entwicklung. Das Produkt von x Faktoren 

hat die Form (1 + d + • • • (1 + 

(1 + <)»= 1 + ^, < + BJ* + 0,«» + ■ • • + HJ'-^ + Kxt'. 

■ * « • ■ 

Der letzte Koeffizient ist gleich Eins. Die Funktion 

Ax hat einen Sinn für a? ^ 1, 

-tßx ff yy ff ff ^ ^ Ä; 

^X ff ff ff ff ^ ^ ^f 



Die Anfangswerte sind 

Nach Multiplikation mit 1 + ^ findet man 
(1 + ^)-+i= 1 + (1 + Äx)t + {Äx+ Bx)t^+ {Bx+ Cxy + 
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Anderseits ist 

(1 + ^)«+i =^ 1 + A^^tt + jB^+ifi + a+i^' + • • • 

Daraus folgt 



Die Auflösung dieser Gleichungen gibt 

Ax= 21 == Ä! + c; es ist c = 0, weil J.^ = 1 ist. 

Bx-^ZAx^Zx^'^^^ + c, 

Aus der Bedingung J?j = 1 folgt c = und 

^ a;(a?-l) 

fi y^T? a;(a?-l)(a;-2) , ^ 
Cx^lJBx ^TgTj + 0. 

Die willkürliche Eonstante c ist = 0, weil C3 == 1 ist. 
In dieser Weise erhalten wir die binomische Entwicklung 

• (i + ty^i + xt + ^-^fi+^^^^^ 

■^ l'2'B...x 

Der letzte Koeffizient ist gleich EinS; wie es auch sein sollte. 

53. Entwicklung von oos xt nach Potenzen von oos f. 

Nach der Moivreschen Formel ist 

(cos t + isrnty^ cos xt + i sin xt 

Der Vergleich der reellen Teile beider Seiten dieser Gleichung 
ergibt 

cos xt = cos*^ — ^^f ~ ^ sin^^ • cos*-^^ 

1 • ä • o • 4 

Wenn man die geraden Potenzen von sin t durch cos t ausdrückt; 
so findet man eine Gleichung von der Form 

t cos xt = Äx cos* t + Bx coB^'^^t + Cx cos*"* t-\ 
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Es sollen die Koeffizienten Äx, B^, 0^, ... bestimmt werden. 
Der Kürze halber wollen wir cos t mit u bezeichnen. In der 
Entwicklung 

cos xt = ÄxU^ + BxU'-^ + Ca,M*-* + D^u^"^ H 

hat die Funktion 

Äat einen Sinn für a? ^ 1, 

Bx n V 3) ^^ ^f 

J^X }} jy w ^ ^ O, 



Aus der Formel cos t^ A^u folgte daß A^ = 1 ist um die 
übrigen Anfangswerte zu bestimmen^ genügt es^ in den Formeln 

cos 2< = ^w*+ B^y 

cos 4/ =^ A^u^ + B^u^ + C4, 

cos 6/ = A^u^ + ^eW* + Ce w^ + A 



/ = — oder w = zu setzen. Man findet 
2 

-Og = 1, 6/4 *= 1 , X/ß =5 1 ; ... 

Aus der Relation 

cos {x + \)t + cos {x — 1)^ = 2 cos t • cos ir< 

folgt die Identität 

Ax^xu^^^ + (:er+i + ^-i)w*-^ + (a+i + Bx-i)w'-» + • • • 

= 2AxU'^+^ + 2J?;,w*-i + 2Ca,w*-* + • • • 
Der Vergleich der Koeffizienten führt zu den Differenzengleichungen 

Bar+l— 2JBa;= — -^x— 1, 
Caj^l — 2 Ca? = — Bx—1, 



Diese Gleichungen vereinfachen sich; wenn man neue Funktionen 
einführt. 

Wenn wir setzen 

A=2*-i.ax, Bx=^-2^-^'lxy Cx^2^-^'Cxy.^.y 

so finden wir ^ 

C^-l-l — C;c = 6«—!, 
<?a?+l— dx=^ Ca:_-i, 
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Die Anfangswerte sind 

Oi^l, 62 = 2, 64=2, ^6=2,... 

Aus der ersten Gleichung folgt; daß ax konstant ist; da aber 

«1=1 ist, so ist 

aaj=l und J.a.==2*"~^ 

Die zweite Gleichung gibt 

6a, = Z 1 =x + c. 

Weil 62 ^^ 2 ist, so ist c = 0, 

ha:=x und JB, = — 2*~^ • ir. 

Die weiteren Bechnungen sind 

c. = Z(a; - 1) = ^ :^ — ^ + c, 

32, . 

(a;-l)(a;-.2) , 

, (a?~2)(a?--3)(a?-A) . .v 

^_^ p-.),.-4)(,-a) _(^_,)j 

__ {x-^)(x-^)(x-6)(x-&) (a? - 5) (a; - 6) 

~" 1.2-3.4 1-2 "*"^' 

<v 54-3-2 32, __ ^ 

"^^1. 2. 3. 4"" 1:2 "^^^ ^""^^ 

_ (j?~3)(a;~4)(a;--6)(a?--6) (a; -- 6) (a? - 6) 
^'"" 1.2- 3.4 1-2 ' 

- ^ '^ r (a?-4)(a;-6)(a;-6)(a?~7) __ (a; ~ 6) (a; - T) ! 

'*"~^L 1-2-3.4 1-2 J 

_ (a;-~4) (a?- 6) (a;~ 6) (a;- 7) (a?-8) (a;- 6) (a;~ 7) (x - 8) 
~" 1-2. 3-4. 6 1-2-3 

c) _ 6-6-43'2 4-32 , _^ 

(a;-4)(aj-5)(a?--6)(a;-7)(a;-8) (a? ~ 6) (a; - 7) (a? ~ 8) 



+ C; 



f^- 



1-2-3-4-6 1-2-3 



Die gefundenen Ausdrücke für Cx, d^, Cx, fx lassen sich ver- 
einfachen. Es ist 
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^-(i-f)(i-|)-i = ^' 



6* = 



123 

(x — b)(x — 6) r/^ *\ /i ^\ il rc(a: — ö)(ir — 6) (a: — 7) 



[(' -!)('- f ) - 



1 • 2 L\ 8/ \ 4/ J 1 • 2 • 3 • 4 

(aj-6)(aj~7)(ir-8) T/^ a;\ /^ x\ ^"1 a;(ir-6)(a:-7)(a:-8)(a:-9) 



f.- ''-T;T'"" [('-i)('-f)-t]° 



12-3 4-5 



f 



Die gesuchten KoefBzienten haben also die Werte 

X=2--S 5,= -2'-».^, 0^:=2--«.^^^, 

T) 9^.-7 a;(a?-4)(a;-6) jp _ ox-9 a?(a? - 5) (g ~ 6) (a? - 7) 

^* * 123 ' ^a:— ^ • 12 -8.4 ' 

TT ^_ oar-ii a?(a;-6)(a;~7)(a?-8)(a?-9) 
* ^ * 12. 3-4. 6 '* 

P. L. Tschebyscheff hat bemerkt, daß die ganze Funktion 
^t«n Grades 

cos X't 

fw == -^i^ir = »* + Pi»'~^ + A«*-* + • • • + i),_iM +p^ 

im Intervall (— 1, + 1) weniger von Null abweicht, als jede 
andre ganze Funktion von derselben Form: 

(Oeuvres de P. L. Tchebychef, t. I, p. 295—301. Vgl. J. Ber- 
trand, Traite de calcul differentiel et de calcul integral, t. I, 
p. 512— 521.) 

Die Abweichung von Null ist der größte absolute Betrag 
einer Funktion. Wenn z. B. 

\F{u)\'^A für a<w<6, 

so ist A die Abweichung von Null der Funktion F{u) im Intervall (a, &). 
Die Abweichung von Null der Funktion 

2 

im Intervall (— 1, + 1) ist gleich ^__^ ' 

Wir wollen voraussetzen, daß im Intervall (—1, +1) die 
Funktion (p(u) weniger von Null abweicht, als f(u). 
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Für 



ist 



^ « 2» {x-'l)7C 

^ X X X 

^ 2ä 2« (rc — 1)ä ^ 

W = 1, cos — ; COS — } ' ' f COS ^ —f — 1, 

' X X X ' 



Wenn wir der Kürze wegen die Bezeichnung einführen 

Mjt= COS 9 (Ä = 0, 1, 2, . . ., x)y 

X 

so bestehen nach nnsrer Yoraussetznng die Ungleichheiten 

fi^o) - VK) > 0, /^(wi) - 9(wi) < 0, /•(%) - 9 (wg) > 0, . . . 
In jedem der Intervalle 

(Mq? Wl); (^1, %), . . ., (Wa;-1, War) 

hat die Gleichung /*(w) — 9)(w) = mindestens eine Wurzel und also 
im Intervall (—1, +1) mindestens % verschiedene Wurzeln. Dies ist 
aber unmöglich, weil die ganze Funktion 

f{u) - 9)(w) = Ol ^ ffi)«**""^ + (i>8 - 38)w"'~^ + • • • + Oa;- 2a:) 

den Grad x—\ nicht übersteigen kann. 

Damit ist also der Satz von Tschebjscheff bewiesen. Diesen 
Beweis verdanke ich einer mündlichen Mitteilung von Herrn 
A. A. Markoff. 

(Vgl. Lacroix, Traite du calcul dififerentiel et du calcul integral. 
Nr. 1052— 1054) 

Drittes EapiteL 
Lineare Differenzengleichnngen mit konstanten Koeffizienten. 

54. Partlkul&re Lösungen der homogenen Olelohong. 

Um die Lösungen der Gleichung mit konstanten Koeffizienten 

ZU finden^ werden wir in die linke Seite dieser Gleichung einsetzen 
(2) y^=a*W:r, 

wo Wa; eine ganze Funktion m*®^ Grades ist. 
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Mit Benatznng der Formeln 

findet man als Besnltat dieser Substitation die Identität 
+ a*+^ [na»-* + (w — l)i)ia'*-* + (» — 2)|>2a»"* H h i>«_i] ^ «*x 






+ a*+«.n(n - l)(n - 2) . . . 2. 1 • jTg-^- 

Zur Abkürzung wollen wir die Bezeichnnng einführen 
(3) f{g) = t^ + p^s!^-^ +PiS^-'^ + ••• + !>„* 

Die gefundene Identitilt nimmt dann die Form an 



(4) 



+ a'+« . f(-)(a) . i4^- 



Unter v ist hier a^'W, zu verstehen. 

Die rechte Seite von (4) verschwindet, wenn a eine Wurzel der 
Gleichung f(js) =» und Ux= 1 ist. Eine partikulare Lösung von (1) 
ist dann a'. 

Wir wollen voraussetzen, daß a eine ^-fäche Wurzel der Gleichung 
f(/s) = ist. Die rechte Seite von (4) ist dann 

««+* . ^(*)(a) . ^^ + a»+*+^ . fi^+^Ka) • fT^^ + ' ' ' 

Dieser Ausdruck ist gleich Null, sobald Ux eine ganze Funktion 
vom (Je — 1)**^ oder niedrigeren Grade ist. Wir setzen 
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Mo = 0, Awa,= l; 



Dann ist 

Uje =1, Wa: = X, Ux == a;(a: — • 1), ... 

Ux = a?(a: — l)(x — 2) . . . (aj •— Ä + 2); 

wir haben also Je partikuläre Lösungen von (1) 

a*, a* • a;, a' • aj(a: — 1), ... , a* • aj(a; — 1) . . . (a; — Ä + 2). 

In Nr. 48 war es schon gezeigt; daß den partikularen Lösungen 
einer linearen homogenen Gleichung ein konstanter Faktor hinzu- 
gefögt werden kann. Deswegen kann man die gefundenen partiknlären 
Lösungen durch die folgenden ersetzen 

a% xa'^'~\ x{x ■— l)a*-^, . . ., x{x — 1) . . . (a? — ft + 2)a*~*"^ 

• 

um diese Lösungen zu bilden^ genügt es, die Funktion ^ zu 
nehmen und Qc - 1) sukzessive Ableitungen nach t zu bilden und 
in den Resultaten t durch a zu ersetzen. 

55. Allgemeine Iiösung der homogenen Oleichong. Bei 

der Auflösung der Gleichung 

(1) y,+n+i>iy«+«-i+Ay.+«-2+---+i>-y.=o 

werden wir zwei Falle unterscheiden: 1. die entsprechende algebraische 
Gleichung 

hat verschiedene Wurzeln und 2. die Gleichung f(g) = hat gleiche 
Wurzeln. 

Wenn wir im ersten Falle die Wurzeln mit 

bezeichnen, dann sind nach Nr. 54 

partikulare Losungen von (1) und nach Nr. 48 genügt 
(5) Vx- ^1«? + C^2«| + C^a'^+'^'+Cna- 

der Gleichung bei willkürlichen Werten der Konstanten 

^v ^fy ^Sf • • •; ^«» 
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um zu beweisen; daß der Ausdruck (5) die allgemeine Lösung 
ist; genügt es zu zeigen, daß 0^, C^, . . .7 Cn sich durch willkürlicli 
gegebene Werte Pq, y^, y^^ . . ., y^^i ausdrücken lassen. 

um das System 

yi^C^a^ +C^ai + Gj«» H h C„a„, 

(6) j y2=C,al +C,al + C^al +... + aaJ, 

^y„«i = C'iaj-H C?2a;-i+ C^a^'''+ • • • + 0„a~-i 

aufzulösen, werden wir die Gleichungen des Systems mit unbestimmten 
Faktoren 3„__i, ^n—v ffn— s^ • • • ? & niultiplizieren und dann addieren. 
Das Resultat nimmt die Form an 



(7) 



wenn man die' Bezeichnung 

einfuhrt. 

Um Cj zu bestimmen, werden wir q^, ?!,... ff„__i so wählen, daß 

9)(ai) = l, 9(0,) = 0, 9(05) = 0,..., q>{an)^0 

wird. Es ist in diesem Falle 



^ (2?-a,)(g~aa)...(g~a^) 



(«1-««) 
und die Relation (7) gibt 

In derselben Weise lassen sich C^, C^, ,,,, Cn durch y^, y^, 
Vif ' • f Vn—i aiisdrücken. 

Damit ist also bewiesen, daß der Ausdruck (5) wirklich all- 
gemeine Lösung von (1) ist. 

Jetzt wollen wir voraussetzen, daß 

gleiche Wurzeln besitzt. 

In allen einzelnen .Fallen ist die Beweisführung dieselbe. Des- 
wegen genügt es, nur einen Fall zu betrachten. Wir werden voraus- 
setzen, daß 



(10) 
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ist und daß die Zahlen »i, «4, a^y ..., a„ voneinander ver- 
schieden sind. 

Da wir n partikuläre Lösungen von (1) kennen (Nr. 54), so 
genügt der Ausdruck 

unsrer Differenzengleichung. Es ist noch zu zeigen, daß das System 

(y„_a = C,aJ-^ + C,(«-l)a;-*+^s(«-l)(«-2)ar' + ^4«;~'+- ' •+^««r^ 
y«-8 = Cf.aj- " + C, («-2) c^-» + (73(«-2)(«-3)«r * + C'*«;-* + • •+ ^X"*' 

yi^Ci«! +Ci-1 +Cg-0 +^40^ +-+C„a„, 

yo = Ci-l +Ci-0 +05-0 +C^■1 +-+C„-1 

in bezug auf C^, Cg, . . ., Cn lösbar ist. 

Wenn man die Gleichungen (10) mit q^, ffi, ..-, 3„_2, Q„_i 
mnltipli2dert und znsammenaddiert, so findet man 

j 2oy„-i + 2iy„-» + • • + ««- j^i + ««-1^0 

Hier ist 

<p(e) = go«"~' + äi^~* + ■ • + 3„_2« + q„_y 

I 

Um Ol zu bestimmen, genügt es, die Koeffizienten Qq, q^, . . ., 
q _2* ^n—i ®^ ^^ wäMen, daß 
(p{a^) = 1, g?'(ai) = 0, 9>"(ai) = 0, (p{a^) = 0, . . . , g?(a„) = 0. 
Um die Gleichungen 

9(^4) = 0, 9)(a5) = 0, . . . , g?(an) = 
zu befriedigen, genügt es, zu setzen 
(12) (p(0) =>(z- a^)(0 — ag) . . . (^ - an)t{d), 

wo il;(z) eine willkürliche Funktion ist. 
Aus (12) und aus den Bedingungen 

folgen die Werte für ^(aj, ^'(a^), ^"(«i). Die Funktion vom 
niedrigsten Grrade, die allen diesen Forderungen genügt, ist 

fp(0) = (;sf - a^{g - Ob) ... (^ - a«) [^(o^) + (^ - «O^'K) 

Seliwanoff, Differenzenrecbnong. 6 
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In derselben Weise lassen sich die Eonstanten C^ und C^ durch 
Vq} tfiy "f y„-.i ausdrücken. 

Um C4 zu bestimmen^ genügt es^ zu setzen 



(p{0) 



(g~fl4)»(g~at)(g-ag)...(g~a^) 
(«4 - öi)"(a4 - «ö) («4 ~ «e) • • • («4- «n) 



Ebenso lassen sich die übrigen Eonstanten C^, C^, ...^ Cn 
bestimmen. 

Damit ist also bewiesen, daß (9) wirklich eine allgemeine 
Lösung der vorgelegten Differenzengleichung ist. 

Durch die Einführung andrer willkürlicher Eonstanten kann 
diese Lösung eine einfachere Form annehmen. Da der Ausdruck 

Ci + C^a^^x + Cj^a-^x(x - 1) 

eine willkürliche Funktion 2*^ Orades ist, so kann man ihn durch 
A^ + A^x + Ä^x^ ersetzen. 

Wir bekommen die folgende allgemeine Lösung unsrer Gleichung 



(13) 



y, = (^ + A« + ^«*)<»? + c^»: + • • • + c,ai- 



Li allen Fällen gibt es also solche n Lösungen u(^>, u^'), . . . , »("' 
der gegebenen homogenen Gleichung, daB man das System 



y,_x = G.u^U + Gt^U + • • • + <^„«Slx 



in bezug auf. C^, Cg, . . ., 
Determinante 



(14) 



<) 
«« 



(7« auflosen kann. Deswegen ist die 



«{*) ... «W 
mW ... «W 



«« 1 »?> 1... Mir', 

von Null verschieden. 

Als Beispiele wollen wir einige Gleichungen mit ihren allgemeinen 
Lösungen angeben. 
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66. G-lelohnngen mit reellen Koeffizienten. Wenn alle 
Koeffizienten der Gleichung 

reell sind^ so kann man verlangen^ daß in der allgemeinen Lösung 
keine imaginären Ghrößen auftreten. 

Wir werden Yoraussetzen^ daß die entsprechende algebraische 
Gleichung 

eine dreifache komplexe Wurzel hat, die wir mit a^ bezeichnen. 
Nach einem Satze der Algebra hat diese Gleichung eine dreifache 
konjugierte Wurzel, die wir mit a^ bezeichnen. Die übrigen Wurzeln 
«,, a„ . . ., o„ soUen reeU tmd voneinander verschieden sein. 
Die allgemeine Lösung nnsrer Gleichung ist 

Wenn 

^1= Q (cos a + i sin a) 
gesetzt wird, so ist 

«4 = (> (cos a — i sin a), 

aj = 9* (cos a;a + i sin xa), 

aj = 9* (cos xa — isin. xa). 
Deswegen ist 

(Ci + G^x + C^x^al + {C^+ C,x + C,x')al 
= [(Ci + C^x + C^x^ + (04 + C,x + C^x'^^Q^ cos ax 

+ [(Ci + ^2^ + C^s^') i -(G^+C^x+ C^x") {\ ^ sin ax. 
Wenn man neue Konstanten einführt 

Ä= C^+C^, -4i= Cj+Oö; -^2= Cq+ Cq, 
( SO erscheint die allgemeine Lösung unsrer Gleichung in der Form 
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y^ = (>* • [(-4 + A^x + A^x^) cos xa + {B + B^x + B^x^) sin xa\ 

+ C.aj + • • • + (^nKf 

die keine imaginären Größen enthält. 
Beispiele: 

y^= 2*- [^ cos If^ + B sin (^x)]* 

y^= Ci cos (ä; J) + (7, sin(a; J) + Gjcos (a;?^^) + C^ sin(Ä;^). 

67. VoUst&ndlge G-lelohnng. Um die Lösung der voll- 
ständigen Gleichung 

(15) y^^n+Piyx+n^i+p^y.+n^.+ '+Pnyz-Q. 

zu suchen^ werden wir voraussetzen, daß die Lösung der entsprechenden 
homogenen Gleichung 

schon bekannt ist. Es sei die allgemeine Lösung dieser Gleichung 

• * 

Statt der willkürlichen Eonstanten werden wir unbekannte 
Funktionen von x einführen, und der Ausdruck 

(16) y, = cwuw + owmW + • ■ • + crw«w 

soll der vollständigen Gleichung genügen. Eine Bedingung haben 
wir zu erfüllen und n Funktionen zur Verfügung. Deswegen können 
wir w — 1 neue Bedingungen einführen. Lagrange wählt dieselben 
so, daß y^ , 1, J/a; 1 2; • . •, J/af-L-n— 1 ^®s®l^® FoHU haben sollen, als ob 
die Cx Konstante wären. 

Da Ge-f-i = Coj + AC« islv so finden wir das System von Gleichungen 

i7 X XX 'XX 'XX' 
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Diese Formeln gelten, wenn die Funktionen Ox den Bedingungen 
genügen 



(17) 



«<1, 

« + 2 



aj-f-w — 1 X x-^n — 1 X x~\-n—\ x 

Es ist 

und es kann nicht mehr y , dieselbe Form haben, als ob die Gx 
konstant wären. 

Wenn wir die Werte y^, y^_^^, y^^^, . . ., y^_^^ in die gegebene 

vollständige Gleichung einsetzen^ so &llen die Glieder mit Cj^^\ 
Cf\ . . . , G^^ weg und es bleibt die Gleichung 

(18) u%^ . ACW + wW „ . AC(«) + . . . + u^Xn ' ^Ci«) « . 

^ '^ aj-f-n a; ' x-\-n x ' ' x-f-n x ^x' 

Aus den n Gleichungen (17) und (18) mit den n Unbekannten A (7^^), 
tikC^\ ..,, A(7^**) lassen sich dieselben bestimmen^ wenn die Determinante 



F{x) = 






^%n »^u • • • «?;« 



von Null verschieden ist. 

Wir haben schon gesehen (Nr. 55)^ daß die partikulären Lösungen 
^x ^ ^x^f • • • ^ ^Ir^ ^^^ homogenen Gleichung so gewählt werden können, 
daß die Determinante 



F(- 1) = 



<) 



mW 



««) 






m(*) . . . w^") 



>t-— 1 « — 1 n — 1 



nicht gleich Null ist. 
Aus der Relation 



u. 



.,+n + Px%+n-^xP2^x+n-^Z+'"+Pn^xK+l+Pn%-^ 



folgt, daß 
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Fix) = 



**«-}- 1 






»^:\.i 
»?;» 






«i?;«-! • • • 



»+«— 1 



"*•*•* -■•«X -■••IX 



weil man zu den Elementen der letzten Zeile die Elemente der vorigen 
Zeüen, resp. mnltipUziert mit p^, j)^, . . ., jp^_„ addieren kann. 
Es besteht also die Relation 

Wäre F(x) = 0, so müßte auch F(x — 1) = sein, weil p^ nicht 
NuU ist. 

Aus denselben Gründen wäre auch 

F(x-2)^0, JF(a?-3) = 0, ..., F{-1) = 0, 

was aber nicht der Fall ist. 

Es ist also. bewiesen, daß F(x) nicht Null ist. Die Gleichungen 
(17) und (18) sihd lösbar. Man findet 

X X ' X X * ' X X * 

Die Summation gibt Resultate von der Form 
Demnach ist die Lösung der vollständigen Gleichung 

X X X X * 'XX 

Es ist die allgemeine Lösung, weil die Determinante (14) von 
Null verschieden ist. 

Diese Methode von Lagrange heißt Variation der will- 
kürlichen Eonstanten. (Lagrange, Oeuvres, t. 4, p. 156. Recherches 
sur les s^ries r^currentes.) 

Als Beispiel nehmen wir die Gleichung 

yx+2-'^yx+i+^yx-^^ 

ACW+6*+^-ACW=0 



Wir setzen 

Es ist 



wenn 

ist. 
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Es wird y^ der Gleicliimg genügen^ wenn außerdem 

X X 

ist. Daraus folgt 

Die Summation gibt 

Die allgemeine Lösung unsrer Gleichung ist also 

68. Man kann noch in einer andern Weise von der Losung 
der homogenen Gleichung zur voUständigen (15) übergehen. Es 
genügt; eine partikuläre Lösung der vollständigen Gleichung zu 
kennen (Nr. 48). Wenn das letzte Glied 

ist, wo V:, eine ganze Funktion ist, dann ist es leicht, eine solche 
Losung zu finden. 

Wir setzen ^ 

wo Ux eine ganze Funktion ist, und führen diese Größe statt y^ in 
die vollständige Gleichung ein. Vermöge der Identität (4) findet 
man die Gleichung 

a'f{a) . Ux + a'+^ • f{a) - Aux + a*+« • f"(a) • -^ + - • • 

.., + a^+n. fin) (^) . ^^"^^ == a' . Vx, 

die zur Bestimmung von Ux dient. 

Der Faktor a* kann gehoben werden und es wird diese 
Gleichung die Form annehmen 

(19) f(a)u^+af'{a)Au,+aY'ia) ■ Yr|+- ' -l-aY<'"(o) • TTäT:^"»- 
Wenn a keine Wurzel der Gleichung 

ist, dann ist u, Ton demselben Grade wie Vx. 
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Es sei 

__ x{x-l){x-2i).,,(x-m + l) x(x-l)...(x-m + 2) 

^x-CCq 1.2.8...W ■*'"l 1.2.3...(w-l) 

, x(x—l)..,(x — m + B) , , x(x-'l) , , 

+ ^ 1.2...(m-2) + • • • + «^-« ~rT~ + "^-^^ + "^^ 

— A a? (a? — 1) (a; — 2) . . . (a; — m -f 1) , ^ a?(a; — 1) ... (a? — m 4- 2) 
^* ~ P^ 1.28.. .w •" ''l 1.2.3...(«-1) 

, a a;(a;-l)...(a;-w + 8) , ^ xjx-l) , ^ ^, /> 

+ Pa 1.2...(m-2) + • • • + Pm-i 1.2 + Pm-i^ + Pm. 

Die Zahlen /3 sind gegeben und die a werden gesncht. 
Aus (19) findet man durch den Vergleich der Koeffizienten 

ccif{d) + a^af(a) = ß^, 
a,f{a) + a,af\a) + a,a^tl^ = ft, 

«srW + «,af («) + «la»^^ + «o^^ff^ » A 

Da f{a) nicht Null ist^ so bestimmen sich sukzessive die Zahlen 
^Q} ^\y ^y ^} • • • 

Wenn a eine %fache Wurzel von f(ßi) = ist, so nimmt die 
Gleichung (19) die Form an 

Die linke Seite dieser Gleichung ist eine ganze Funktion 
m*®^ GradiBS, wenn Ux vom (Ä + m)^^ Grade ist. 
Man kann setzen 

a? (a? — 1) . . . (a? — Ä; — w + 1) aj (a; — 1) . . . (a? — ifc — w + 2) , 

^^ "" ^ 1.2.3...(Ä; + w) ^ ^1 1.2.3...(Ä; + w-l) ' 

, a5(a;— 1). . . (a; — Ä;) a;(a!; — 1) . . .(a; — A;-f 1) 

"^ ^-^ 1.2.3... (Ä+l) "*" ^"* 1.2.3...Ä; ' 

und zur Bestimmung der Koeffizienten findet man die Gleichungen 

**o" 1.2...&~P0' 

"1*1.2...* + '^" 1-2.8...(Ä! + 1) Pi' 

**«" 1.2...*+'^* 1.2...(fc+l) + *^* 1.2...(ifc + 2) P«- 
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Diese Gleichnngen sind lösbar^ weil a und f^^\d) nicht gleich 
Null sind. 

Z. B. wollen wir eine partikuläre Lösung der Gleichung 

suchen. Hier ist a = 1 und dieser Wert ist eine einfache Wurzel 

der Gleichung 

^2-7^ + 6 = 0. 

Deswegen gibt es eine partikuläre Lösung von der Form 

V — xjx — l) 

Jlx — CCq j-y-r [- CC^X, 

In die gegebene Gleichung setzen wir ein 

Yx= bx + cx^ 

und wir finden 

I. 3 1 

* = 6Ö' ^ = -10- 

Nach dem bekannten Satze (Nr. 48) ist die allgemeine Lösung 
dieser Gleichung ' 

Die Gleichung ^ 

hat eine partikuläre Lösung von der Form 
Aus der Gleichung 

folgt; daß cCq^ — ist. Deswegen ist die allgemeine Lösung unsrer 
Gleichung 

Wir wollen jetzt einige Anwendungen der Differenzengleichungen 
höherer Ordnung zeigen. 

59. Bekurrente Beihen. Dies sind solche Reihen 

Vof Vi) Vtf ' ' *f Vxy Vx-^v ' • •' 
für die je n aufeinander folgende Glieder einer Relation 

yx+n + Pxyx:{.n-^l + P2yx+n^2+-^+Pnyx-^ 

mit konstanten KoefOzienten genügen. 

Seliwanoff, Differenzenrecbnung. Q* 



« « > 
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Die Auflösung dieser Differenzengleichung dient dazu^ um das 
allgemeine Glied y^ als Funktion von x darzustellen. 
Als Beispiel wollen wir die Reihe betrachten 

0, 1, 1, 2, 3, 5, 8,..., 

deren jedes Glied die Summe der beiden vorhergehenden Glieder ist. Je 
drei nacheinander folgende Glieder sind durch die Relation verbunden 

Die Anfangswerte sind 

Die Wurzeln der. quadratischen Gleichung 
sind _. ,-. 

Deswegen' ist die allgemeine Lösung der Differenzengleichung 

Wenn C^ and (7, durch die An&ngswerte y, und y^ bestimmt 
werden, findet man, daß das allgemeine ÖUed unsrer Reihe 

^.=i^[(^r-(^)'] 

ist. 

Die Potenzreihe 

F(t) = j/o + Vi^ + %<' + j/s^" + • • + yj' + • •' 

deren Koeffizienten eine rekurrente Reihe bilden^ ist gleich einer 
rationalen Funktion von U Die Relation zwischen den n Koeffizienten 
der Reihe sei 

Nach Multiplikation von F(() mit 
findet man 

yo + (s/i+Ayo)^ + (y« +Piyi+p%y^)^^ + {vz ^PiVt +p%yi +Ayo)<* + 
• • • + (y„«i +jPiy,^, +i>2y«-8 + • • • +i>„-iyo)*~"'- 

Die Koeffizienten der folgenden Potenzen sind die Werte der 
Funktion 

yx-^n+Piy.+n--i+p^y.+n^t+'"+Pny. 

für rr = 0^ 1^ 2^ . . . Deswegen sind alle diese Koeffizienten gleich Null. 



J 
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Umgekehrt^ die rationale Funktion 

laßt sich in die Reihe 



• • • 



entwickeln, für die je n aufeinander folgende Koeffizienten durch die 
Relation 

verbunden sind. Die n ersten Koeffizienten bestimmen sich aus den 
Gleichungen 



Also können die Differenzengleichnngen dazu dienen, um die 
Koeffizienten der Entwicklungen rationaler Fimktionen in Potenz- 
reihen zu bestimmen. 



60. TBOhebysoheSteher Beweis eines Satzes von Lain6. 

In seinen Vorlesungen zeigte Tschebyscheff die folgende An- 
wendung der Differenzengleichungen. 

Man soll bestimmen^ wie groß die Anzahl der sukzessiven 
Divisionen sein kann, um den größten gemeinschaftlichen Teiler 
zweier Zahlen a und b zu fiinden. 

Es sei 

a=^bqi + r^, b-=r^q^ + r^, . . . , r„-_2 = ^n-i?„ + ^«, r« == 0. 

Die Anzahl der Divisionen ist hier n. 

Weil alle Quotienten größer oder gleich Eins sind, so ist 

a^b + r^, b^r^ + r^, r^^r^ + r^, . . ,, 

Daraus folgen die Beziehungen 

Tn^O rn-i^l, r„_2^1, rn^i^2f 
rn^A ^ 3, r„_5 ^ 5, rn-e ^8, ... 
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Wenn wir die Zahlen 

^0 = 0, yi==i, % = !; y;« = 2, ..., y^=y,_i + y,_2, ... 

einführen^ so ist 

Wenn für einp. ganze positive Zahl m die Ungleichheit 

besteht, dann ist notwendig y^ > y^ und m> n. Die gesuchte 
Anzahl der Divisionen ist also kleiner als eine solche Zahl m, für die 

ist. Den Ausdruck für y^ kennen wir schon aus Nr. 59. Es ist 
deswegen die Ungleichheit 

aufzulösen. 
Da 

/5 = 2,2...; 1:0:^,== 1,6...; ^ 0,6..., 

SO können wir y aniuUiemd in der Weise ausdrücken 



i_(i+y± 



y± 0, 



WO 6 zwischen und 1 ist. 
Wenn 

ist; dann ist 

y„^h + (l±e) und folgUch y^>b. 

Es genügt 

log [2, 2 .(6 + 1)] 
= log 1, 6 

ZU nehmen. Wenn die Zahl b aus k Ziffern besteht, dann ist der 
gemeine Logarithmus von 2, 2.(6 + 1) ungefähr gleich Je, log 1, 6 

ist beinahe — und 

Daraus folgt der Satz' von 0. Lame: Die Anzahl der Divi- 
sioneU; die auszuführen ist, um den größten gemeinschaft- 
lichen Teiler der Zahlen a und 6 (a> 6) zu finden^ ist kleiner 
als die fünffache Anzahl der Ziffern der Zahl b. 

Lame hat diesen Satz auf elementarem Wege bewiesen. (Paris. 
Comptes rendus, 19. [1844], p. 867.) 



Dmck von B. G. Teubner in Dresden. 



r 



Verlag toh B. Gr. Tenbner ia Leipzig. 

Encyklopädie der 

Mathematischen Wissenschaften 

mit Einscbliifi ihrer Anwendungen. 

Herausgegeben im Auftrage 

der Akademien der Wissenschaften zu München und Wien und der 

GeseUschafl der Wissenschaften zu ööttingen, 

sowie unter Mitwirkung zahkeicher Fachgenossen. 

In 7 Bänden, gr. 8. geh. 

Band I: Arithmetik und Algebra, redigiert von W» Fr. Meyer in Königsberg. 

— 11: Analysig in 2 Teilen H« Bnrkhardt in Zürich. 

— III: Geometrie in 3 Teüen Yf* Fr. Meyer in Königsberg. 

— rV: Mechanik in 2 Teüen F. Klein in Göttingen. 

V: Physik in 2 Teilen A. Sommerfeld in Aachen. 

— VI, 1: Geodäsie und Geophysik . . . . E. Wiechert in Göttingen. 

— VI , 2 : Astronomie K. Schwarzschild in Göttingen. 

— YII: Historische, philosophische und didaktische Fragen behandelnd, sowie 

Generalregister. (In Vorbereitung.) 

Bisher erschien: 

I. Band. 1. Heft. 1898. n «J 3.40. 2. Heft. 1899. n. Jl S.40. 8. Heft. 1899. n. «J 8.80. 4. Heft. 1899. 
n. Ul 4.80. 5. Heft. 1900. n. Ul 6.40. 6. Heft. 1901. n. «J 7.80. 7. Heft. 190:^. n. Ul 8.60. 
II. — I. Teil. 1. Heft. 1899. n. Jl 4.80. 2/8. Heft. ISOO. n. Ul 7.60. 4. Heft. 1900. n. Ul 4.80. II. Teil. 
1. Heft 1901. n. «J 5.80. 
m. — IL Teil. 1. Heft. 1908. n. ^4.80. IQ. Teil. 1. Heft. 1903. n:ur5.40. 8/8. Heft. 1908. n. ^16.80. 
IV. — L TeU. 1. Heft. 1901. n. «J 3.40. 8. Heft. •1903. n. «J 4.60. 8. Heft. 1908. n. «J 4.60. U. TeU. 
1. Heft. 1901. n. «J 8.80. 8. Heft. 1908. n. Ul 8.80. , 

V. — I. Teil 1. Heft. 1908. n. U» 4.80. 
TTnter der Presse * 
n. Band. I. Teü. 6. Heft. — V. Band. n. TeU. 1. Heft. — VT. Band. L Teü. 1. Heft. 

Encyidopädie der Elementar-Maüiematilc. 

Ein HandbQcb fOr Lehrer und Studierende. 

Von 

Heinrich Weber, und Josef Wellstein, 

Professor in Strasburg Professor in Gießen. 

In 3 Bänden. Band I (XTV n. 446 S.) gr. 8. 1903. In Leinw. geb. n. Jf 8.— 

Das Werk, dessen erster Band soeben erschienen ist, richtet sich in erster 
Linie an die Lehrer, die darin Anregung finden sollen, ihren Unterrichtsstoff 
auszuwählen und namentlich in den höheren Erlassen zu vertiefen, sodann aber 
auch an Studierende, die eine Anlehnung an die Elemente und Auffrischung und 
Ergänzung früher erworbener Kenntnisse suchen. 

Durch das Zusammenwirken mehrerer Grelehrter hoffen die Herausgeber, die 
möglichste Vollständigkeit zu erreichen. Der erste Band umfaßt den algebraisch- 
analytischen Teil. Der zweite Band, der Ostern 1904 erscheint, wird die Geometrie 
nach ihren verschiedenen Seiten behandeln. Ein dritter Teil, der im Herbst 1904 
zur Ausgabe gelangen soll, wird die Anwendungen bringen, deren St^ff aus der 
darstellenden Greometrie, der Mechanik, Physik und Wahrscheinlichkeitsrechnung 
entnommen ist. Die Vorarbeiten sind so weit gediehen, daß die Vollendung des 
ganzen Werkes im nächsten Jahre erwartet werden darf. 



